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Chapter 1 


Introduction 


This book is designed to help you understand turbomachinery. It aims to help you over some of the 
difficult initial concepts so that your work or study with turbomachinery will be much more fruitful. 
It does not tell you how to design a turbomachine but instead aims to make your other studies, lectures 
and textbooks which go into more depth make much more sense. For those readers not concerned with 
turbomachinery design it might provide all the background they need. It is based on an introductory 
course taught at Durham University for some years. 


There are actually only three really difficult ideas in this book: understanding the cascade view 
(Chapter 1), velocity triangles (Chapter 2)and the concept of reaction (Chapter 6). Once you have 
mastered those three concepts Turbomachinery actually becomes relatively straightforward! 
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This book is available on-line and any comments or suggestions about the book are gratefully 
received by the author. He can be contacted via e-mail at: g./.ingram@durham.ac.uk 


1.1 How this book will help you 


The book is designed to provide guidance on the basics. So if someone is presenting a velocity 
triangle which you do not understand or you have absolutely no idea what a stator is this book will 
help. Armed with this understanding you can then go on use the more complex texts effectively. 


The book also contains examples which illustrate how understanding these basic concepts lead to 
an immediate appreciation of why machines look the way they do. So for example you will rapidly be 
able to see why wind turbine blades are twisted, why a the blade height in a steam turbine increases 
towards the low pressure end and why pump blades often point away from the direction of rotation. 


The most valuable learning experience however is to actually manipulate the ideas contained in 
this text. A series of problems are provided at the end of each chapter with numerical answers - to 
fully understand the material in this book you should attempt these problems. 


1.2 Things you should already know 


This book is directed at readers with a basic knowledge of Fluid Mechanics and Thermodynamics. 
In order to make best use of the book you should have some knowledge of the steady flow energy 
equation, static and stagnation conditions, the perfect gas law, how to use steam tables and charts and 
an understanding of the boundary layer. 


1.3 What is a Turbomachine? 


A turbomachine is a device that exchanges energy with a fluid using continuously flowing fluid and 
rotating blades. Examples of these devices include aircraft engines and wind turbines. 


If the device extracts energy from the fluid it is generally called a turbine. If the device delivers 
energy to the fluid it is called a compressor, fan, blower or pump depending on the fluid used and the 
magnitude of the change in pressure that results. Turbomachinery is the generic name for all these 
machines. 


Somewhat confusingly the word turbine is sometimes applied to a complete engine system on an 
aircraft or in a power station, e.g. “a Boeing 747 is equipped with four gas turbines for thrust”. A 
glossary is in Appendix A on page 137 at the end of the book to help you navigate your way through 
the turbomachinery jargon. 


Turbomachinery is essential to the operation of the modern world. Turbines are used in all sig- 
nificant electricity production throughout the world in steam turbine power plants, gas turbine power 
plants, hydro-electric power plant and wind turbines. Pumps are used to transport water around mu- 
nicipal water systems and in homes, pumps and turbines are also essential in the transportation of 
fuel oil and gas around pipe networks. Gas turbine engines are used to power all large passenger 
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Figure 1.1: Applications of Turbomachinery 


aircraft either in the form of turbo-prop or turbo-fan engines and also through a gearbox they power 
all helicopter engines. 


In short turbomachinery is all around you and is an area worthy of further study! Figure 1.1 shows 
four important applications of turbomachinery, in the top left gas turbine propulsion for aeroplanes, 
in the top right wind turbine power of electricity production, in the bottom left the rotor of a steam 
turbine for power production and a water pump is shown in the bottom right. 


1.4 A Simple Turbine 


There are many variants of turbine, here we describe the operation of a simple turbine so you get a 
feel for what is going on. An outline of a turbine is shown in Figure 1.2. From this view all we know 
about the device is that flow goes into it and as if by magic the shaft turns and produces a torque. 


If we look at the device in an exploded view (Figure 1.3) we see that as well as a number of covers 
and bearings there is a row of aerodynamically shaped objects that don’t move followed by a row of 
aerodynamically shaped objects that provide the torque to the shaft. 


The objects are known various as blades, buckets, nozzles, aerofoils or airfoils. In this book we 
will generally refer to them as blades. The row of stationary blades is known as a stator and the row 
of rotating blades connected to the output shaft is known as the rotor. 


The basic mechanism of operation is as follows (Figure 1.4): 


1. the fluid flows directly into the device in an axial direction (in line with the machine) 
2. the stator blades turn the flow so that it is lined up with the turbine blades 


3. the turbine blades turn the flow back towards the axial direction and turn the output shaft. 
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The key point is that the power extraction from the fluid arises from turning the flow. More 
complex turbines use more than one row of rotors and stators, but all work on the same essential 
principle. A question often asked at this point is that since all the power comes from the rotor can you 
do without a stator? The answer is yes! Wind turbines extract power from fluid with the need for a 
stator. However for flows with much larger energy densities such as those in aircraft engines adding 
a stator allows you to get much more energy out of the subsequent stator row - the reason for this is 
found in Chapter 5. 


1.5 The Cascade View 


There are two key views of turbomachinery used throughout this book (and in turbomachinery design 
in general). These are the cascade view and the meridional view. 


The cascade view arises from looking at the stator and rotor of the simple turbine shown earlier 
(Top half of Figure 1.5) if you look closely at the topmost part of the turbine you can see the blades 
of the stator and rotor outlined in plan view. This is highlighted by a red box. You can actually do 
this for any rotor/stator blade combination around the circumference of the turbine. The fact that 
you can do this for every blade suggests that the plan view may be an excellent way of analysing the 
performance of the machine. 


The 2D cascade view of the simple turbine is shown in the lower half of Figure 1.5. The cascade 
view with a single stator and rotor blade is highlighted with a red box. The relation between the 2D 
cascade view and the 3D real turbine should be obvious. The rest of the cascade view is made up of 
plan views of the other stator and rotor blade combinations. When looking directly down onto the red 
box in the the 3D view of the turbine the movement of the rotor blade appears to be simply from left 
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Rotation 


Figure 1.2: A Simple Turbine 
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Figure 1.3: A Simple Turbine: Exploded View 
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Figure 1.4: Simple Turbine Operation 
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Figure 1.5: Cascade View 
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Very large radius | 


Figure 1.6: The Cascade View as a Large Radius Machine 


to right. So in the cascade view the rotary motion in the 3D model becomes 2D linear motion in the 
cascade view. 


We can then analyse how the turbine blades influence the flow by looking at this 2D cascade view, 
since the cascade view is the same for every blade passage around the circumference of the turbine. 
Although we have completed this for the top of the turbine we can repeat the exercise at any radius 
from the hub of the machine to the tip. 


An alternative way of looking at the cascade view is to say that we are examining an infinite 
radius machine. Consider Figure 1.6 which contains three views, the first is a 3D view of a simple 
turbine, the second shows a sketch of the turbine as viewed from upstream with the blades and hub 
shown in schematic form. To form the cascade view we can approximate the real turbine rotating at 
speed w with a tip radius R = 0.15m and a spacing between the blades of s with a machine with an 
infinite radius and the same blade spacing (or pitch) of s. The rotation of the machine w is replaced 
by a linear motion of magnitude wR where R is the radius of the original machine. 


The actual cascade view involves looking down from the casing to the hub so you get a plan view 
of the blades. Note that in the real machine the pitch s gets larger with larger radius r so the cascade 
view only accurately represents the machine at a single radius. For machines with very large changes 
of radius such as wind turbines we can draw a number of cascade views at different radii. 


The cascade has two “analysis stations” associated with it at inlet and outlet. A consequence of 
the cascade view is the properties of the fluid (pressure, temperature etc) going through the machine 
are assumed constant in the tangential direction since there is no change in geometry or flow between 
one blade and the next in that direction. In the real machine this assumption represents properties 
being constant around the circumference of the machine so that a single value describes the fluid state 
around the whole machine. Analysis stations can also be applied to parts of the turbomachine that 
don’t always have rotational symmetry such as the inlet or the exit pipe - what is assumed there is 
that a single value accurately represents the flow in the inlet or the exit. 
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1.6 The Meridional View 


The meridional view is much more straightforward than the cascade view and is illustrated in Figure 
1.7. On the left of Figure 1.7 is the familiar 3D view of our simple turbine. For the meridional view 
instead of looking at the tip of the blade this time we take a side on view of the whole turbine and 
look at a cross section of the machine at the hub and tip radius. This is highlighted by a red box. On 
the right of Figure 1.7 is the actual meridional view which shows the stator followed by the rotor in 
cross section. The actual machine radius r is usually very large compared to the blade height b and 
so the axis of rotation is not always shown in the meridional view. 


1.7 Assumptions used in the book 


It is easy to see how the real turbomachinery flow field is three dimensional and unsteady now that 
the complex geometry of machine has been shown. In addition the flow is compressible so density 
changes have to be accounted for. However to introduce the basic concepts we can dispense with a 
great deal of this complexity by making a number of assumptions about the flow field. 


1. The flow is symmetric in the circumferential direction. There is no variation in the flow from 
one side of the blades to another 


2. We consider a mean flow (technically called a stream surface) between the hub and casing. 
This is reasonable for short blades, for longer blades the “trick” is to repeat the calculation at a 
number of radii. 
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Figure 1.7: Meridional View 


3. Flow is steady. Although state of the art blade design requires a consideration of unsteady flow 
most of the turbomachinery in use today has been designed with this steady flow assumption. 


4. Flow follows the blade exactly. There is no deviation between the direction that the blades are 
pointing and the direction that the fluid travels in. (In turbomachinery jargon: the flow follows 


the metal angle of the blades) 


These assumptions may seem quite limiting but most of them are used in the preliminary design 
of all turbomachinery in use today so actually get you a surprisingly long way! 


1.8 Problems 


1. Explain why a bicycle pump is not classified as a turbomachine. 


2. Sketch the cascade and meridional views for a horizontal axis wind turbine such as the one in 


the top right of Figure 1.1. 
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Chapter 2 


Relative and Absolute Motion 


One of the key concepts in turbomachinery is understanding how the flow appears from the point of 
view of components that are rotating compared to those that are stationary. Once this is understood 
this the shape of turbomachinery becomes much easier to understand! Viewing flows from the point 
of view of a rotating component is known as being in the relative frame of reference and viewing flows 
from the point of view of a stationary observer is called being in the absolute frame of reference. 


We start therefore with a simple explanation of relative and absolute motion before ending this 
Chapter with a discussion of how this relates to turbomachines. 
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V 
—— 
W=V-U 


Figure 2.1: Relative and Absolute Velocities for a Cyclist 
2.1 1D Motion 


Consider the everyday activity of riding a bicycle with three cases one where there is no wind, the 
second with a tail wind and a third with a head wind. This is shown in Figure 2.1. The velocity of the 
bicycle we shall label U and call it the “frame velocity”, the velocity of the wind we label V and call 
this the “absolute velocity”. Clearly the absolute velocity V is the velocity that will be experienced by 
an observer watching the cyclist. The wind velocity experienced by the cyclist is called the “relative 
velocity” and given the symbol W. 


The first case shown at the top of Figure 2.1 shows the simplest case, if there is no wind the ob- 
server watching the cyclist will experience no wind and the cyclist will experience a relative velocity 
that is equal and opposite to that of the speed at which he or she is cycling. So the relative velocity 
W=-U. 


The second case concerns a tail wind that is roughly equal in magnitude to the speed of the bicycle 
U. This is shown in the middle of Figure 2.1. In this case a stationary observer would experience the 
wind velocity but since the cyclist is moving at the same speed as the air the relative velocity W will 
be around zero and the cyclist will experience no wind. 


The third case concerns a head wind that is again roughly equal to the velocity U of the bicycle 
in magnitude but not in direction. This is shown at the bottom of Figure 2.1. A stationary observer 
would experience the same wind velocity as in the second case but in a different direction. The cyclist 
however has a very different experience. The relative velocity is made up of their own speed — at (that 
of the first case) added to that of the oncoming wind y. By inspection we can see that W=V-U. 
Since V is negative the cyclist now has to work much harder to maintain the same forward speed. 


This suggests a generalisation of the relationship between relative and absolute velocity: 


V=U+W (2.1) 
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V=0 
—— ee 
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Figure 2.2: Velocity Triangles for an Aircraft Landing 


Or in words absolute velocity is the vector sum of the frame velocity and the relative velocity. A 
trivial rearrangement returns us to the relationship seen in Figure 2.1. 


W=U =V (2.2) 


2.2 2D Motion 


We will apply our new found key rule (Equation 2.1) to one other non-turbomachinery situation to 
illustrate how it works. This situation is one where the motion is in two dimensions. Consider the 
plan view of a aircraft and a runway in Figure 2.2. In the first situation (top of Figure 2.2) there is no 
atmospheric wind V = 0 and so the aircraft simple lines up with the runway and lands. 


The second situation (lower part of Figure 2.2) is where there is a substantial cross-wind, in this 
case imagine that the wind is entirely perpendicular to the runway. What relative velocity (W) does 
the aircraft have to fly at to ensure that the movement of the aircraft (the frame velocity U) results in 
the aircraft arriving on the centre-line of the runway? 


The frame velocity we know is given by the desired path of the aircraft, that is directly towards the 
runway and the absolute velocity is given the atmospheric conditions. The relative velocity is given 
mathematically by the application of our key rule, Equation 2.2. But what if we wanted to sketch out 
the vector? This enables us to understand the direction the aeroplane should be facing. 


To do this we need to use a tool known as a velocity triangle one of the fundamental tools of 
turbomachinery analysis. First we review some very basic vector addition and subtraction rules, 
shown in Figure 2.3. 


e To add two vectors A + B graphically: place them nose to tail and the result is given by 
movement from the tail of the first to the nose of the second. 
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Figure 2.3: Graphical Addition and Subtraction of Vectors 


e To subtract two vectors A— B graphically: reverse the direction of B then proceed with addition 
of vectors as before. 


To apply this to the example of our aircraft we apply the key rule and our knowledge of how to 
put vectors together to end up with the required relative velocity. This is shown in the lower portion 
of Figure 2.2, first the frame velocity U is reversed in direction to form —U, this is then added to V 
by putting them nose to tail. We then draw the line between the start of the vector V and the end of 
the vector —U which gives the relative velocity W. 


This explains why aircraft landing in cross-winds often have to approach the runway at an an- 
gle. If you have an active web connection there are some spectacular examples of this on YouTube: 
http: //uk.youtube.com/watch?v=GHrLB mlir4 


Note that we formed the relative vector W by drawing V then —U but we would end up with the 
same result if we drew the triangle with —U then V. 


All this may seem obvious but it is vitally important before we move onto turbomachinery that 
you are confident in how to draw a 2D vector and how to add and subtract vectors graphically. 


2.3 Velocity Triangles in Turbomachinery 


In this book we consider a Cartesian coordinate system consisting of an axial x, radial r and tangential 
O set of coordinates. The velocity of the frame of motion is denoted by U, velocities in the frame of 
motion are denoted with W and absolute velocities are denoted with V. Consider a turbine consisting 
of a stator and a rotor, the cascade and meridional views are shown in Figure 2.4 along with the 
coordinate system. 


There are three points that are of interest to us entry to the stator, the gap between the stator 
and the rotor and exit from the rotor, these are labelled 1,2 and 3 respectively in Figure 2.4. The 
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Figure 2.4: Cascade and Meridional Views of a Turbine Stage 


combination of rotor and stator is called a “stage” in turbomachinery jargon. These points are the 
analysis stations referred to in Chapter 1. 


At point 1 we have an incoming velocity but as the stator is not moving there is no relative motion 
between the incoming flow and the stator so there is no velocity triangle to draw at this point. 
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At point 2 the flow leaves the stator and enters the rotor. Here there are two frames of reference, 
the flow viewed from the point of view of the stator and the point of view from the moving rotor. A 
velocity triangle can be drawn here. The rotor in the cascade view is moving with a linear (tangential) 
velocity of magnitude wrm where w is the rotational speed of the machine and rm is the mean radius 
of the blades. 


At point 3 the flow leaves the rotor and exits the stage. Again there are two frames of reference, or 
points of view for the flow. That found by viewing from the moving rotor and that found by viewing 
from outside the rotor where there is no motion. 


We can now draw the velocity triangles for point 2 and point 3 in the stage, this is shown in Figure 
2.5. The methodology for this is as follows: 


1. Draw the flow that you know 
2. Draw the blade speed 
3. Close the triangle with the remaining vector 


4. Check that the key rule applies: V=U+W 


This methodology is important and we will refer to elsewhere as the “four step rule”. 


So for station 2 the flow that we know is the absolute velocity at exit from the stator, V. This is 
the flow that we know, recall from Chapter | that the flow follows the metal angle of the blades so if 
with a sketch of the stator the absolute velocity may be drawn directly. To get the velocity triangle 
draw the absolute velocity vector V, draw the blade speed U and then close the triangle with the 
relative velocity W. The result is in Figure 2.5. The final (and vital!) step is to check that the correct 
triangle has been obtained by following the blade speed and relative velocity vectors. If we end up in 
the same place as if we had followed the absolute velocity vector the triangle is correct. 


For station 3 the flow that we know is the relative velocity at exit from the rotor, wW again this is 
because the flow follows precisely the path of the blades and since the rotor blades are moving the 
flow that we sketch on the rotor must be the relative and not the absolute velocity. Having drawn W 
we draw the blade speed U and then close the triangle this time with the absolute velocity V. Again 
the final vital step is to check that the we have the correct triangle by making sure the key rule applies. 
The correct velocity triangle is shown at the bottom right of Figure 2.5. 


Once we have drawn the velocity triangles we can then carry out a series of calculations on 
the fluid going through the turbomachine and the end result might indicate that our sketch is not 
entirely accurate, i.e. the blade speed is much greater than what we have drawn in Figure 2.5 - this 
does not matter! If accurate velocity triangles are required they can always be drawn again once the 
calculations are complete. 


2.4 Velocity Components 


For each velocity in the cascade view we can decompose into axial and tangential components and 
can also express each vector as a magnitude and direction. Axial components are denoted with the 
subscript x and tangential components are denoted with subscript 0. Angles can be measured in a 
number of directions but in this book the axial direction is chosen. The angle made by the absolute 
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Figure 2.5: Velocity Triangles for a Turbine Stage 


velocity with the axial direction is called a and the angle the relative velocity makes with the axial 
direction is called 8. Angles are positive in the direction of rotation. Therefore velocities can be 
specified as a vector V or a magnitude and angle, V and a. 


The one complication to this is that in this book we also deal with machines where the flow has 
a significant radial component - in that case we draw velocity triangles in the radial (subscript r) and 
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Figure 2.6: Velocity Triangles at Station 3 of a Turbine 


tangential plane and angles are measured from the radial direction. This will become clearer when 
radial and centrifugal machines are explained in Chapter 4. 


The velocity triangle at station 3 in Figure 2.5 is shown in Figure 2.6 with the various components 
labelled, in order to indicate that we are dealing with station 3 a subscript 3 is added to all the symbols. 
The relative and absolute flow angles a and 8 are also shown. 


From basic trigonometry the follow relationships apply for any station in a turbomachine. 


vV? = V +V (2.3) 
Vz = V cosa (2.4) 
Ve = Vsina (2.5) 
Vo = V, tana (2.6) 
wW? = Wi+Wwe (2.7) 
W, = W cosß (2.8) 
We = Wsin (2.9) 
Wọ = W,tanG (2.10) 


Aside from trigonometry we can also work out that Wy = Vz for all turbomachinery. The reason 
for that is if we look the basic geometry of a turbomachine such as that shown in Figure 1.4 we see 
that there is no motion of the machine components in the axial direction. That is the stator and the 
rotor remain the same distance apart when the machine is operating. The only time we would have 
movement between the stator and the rotor is if the device had suffered some sort of catastrophic 
failure - there is no normal operating procedure where the gap between the rotor and stator would 
increase! 


Example Consider an Office Desk Fan. It rotates at 200 rpm and has a diameter of 30 cm. Air 
enters the fan at 3 m/s, parallel to the axis of rotation. Calculate the relative velocity (W) at the tip 
of the fan. 
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Figure 2.7: Velocity Triangles for a Desk Fan 


Solution The information given is the absolute velocity V which is the air entering the fan. The 
plan is to work out the frame velocity and thus determine the relative velocity. 


To do this sketch the fan and velocity triangle in Figure 2.7. A sketch of the desk fan is on the 
left hand side of Figure 2.7 and the cascade view and the velocity triangles are shown on the right. 
It is very good practise to sketch the object that you are trying to do calculations on to determine the 
overall layouts - even if (as in Figure 2.7!) the sketch is only approximate. 


The frame velocity U we can obtain from the rotational speed and the radius: 


2m 0.3 
U = wr = 200 x 60 * 2 3.14 m/s 


Trigonometry gives the magnitude of the relative velocity: 


W = VV? + W2 = 33 43.142 = 4.34 m/s 


and the angle: 


V tan 8 = —U B = tan”! (=) = tan”! (=*) = —46.3° 


p is negative as the angle is opposite to the direction of rotation. Many students write V tan G = U 
which gets the correct magnitude but a careful inspection of Figure 2.7 will reveal that this has the 
wrong sign. 


2.5 Problems 


1. An aeroplane approaches a runway at 77 m/s with a crosswind of 15 m/s. What angle does 
the aeroplane have to face into the wind to travel directly towards the runway? Answer: 11° 
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2. An office desk fan rotates at 200 rpm. Air enters the fan at 3 m/s, parallel to the axis of 
rotation. Calculate the relative velocity (W) at the hub of the fan if the hub diameter is 10 cm. 
Answer: 3.18 m/s, —19.2° 


3. The flow at exit from a turbine stator row has a velocity of 100 m/s at an angle (a2) of 70° 
to the axial direction. Calculate the tangential and axial velocity components. The rotor row 
is moving with a velocity of 50 m/s. Calculate the velocity magnitude relative to the rotor 
blades at inlet and the relative inlet flow angle (G2). At exit from the rotor row the relative 
flow angle (83) is —60°. Assuming that the axial velocity is constant across the row, what is 
the absolute exit velocity magnitude and direction? Answer: 94.0 m/s, 34.2 m/s; 55.7 m/s, 
52.1°; 35.4 m/s, —15.1° 


4. For the turbine above, assuming that the relative flow at exit from the rotor row is unchanged, 
calculate the blade speed that would give absolute axial flow at exit (i.e. no swirl) Answer: 
59.2 m/s 
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Chapter 3 


Simple Analysis of Wind Turbines 


This chapter provides an immediate application of the principles of relative motion by using the 
horizontal axis wind turbine as an example. Such a wind turbine is shown in Figure 3.1 as can 
be see from Figure 3.1 the blades are far apart so the influence between them is very small. In 
turbomachinery jargon the pitch s is very large. The interactions between the different blades can be 
ignored in a simple treatment. The wind turbine is one of the only examples in turbomachinery where 
each blade can be considered in isolation and this is one of the reasons that a simple analysis is easy. 


Consider one of the turbine blades shown in Figure 3.1, the speed of each of the three blades will 
be the same and assuming that the wind does not vary over the area of the machine, an analysis need 
only be conducted on a single blade and multiplied as necessary. The interaction of the blades and 
the tower is ignored but since this occupies a small fraction of the 360° rotation of the turbine this 
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Figure 3.1: Wind Turbine Picture and Sketch 


will not influence the analysis greatly. 


Consider what happens if an observer was positioned on the turbine blade about half way along 
the span and a virtual cut made through the blade. If the observer looked towards the hub of the 
machine, they would see an aerofoil profile rotating around the hub. Since the radius of the machine 
is large the rotational motion of the turbine blade can be approximated by a linear motion - in much 
the same way as the earth is round but looks flat as the radius is very large. Such a view of a turbine 
blade is shown in Figure 3.2. The rotational motion of speed w is translated to a linear motion of wr 
the tangential velocity of the rotating blade. If the turbine is facing into the wind the incoming wind 
V will be perpendicular to the rotating blade. 


The velocity triangle for the wind turbine blade can then be drawn according to the four step 
procedure and is shown on the left hand side of Figure 3.2. The flow that is known is the incoming 
wind velocity V which is in the absolute frame of reference. The blade speed, U is then drawn and 
the triangle is closed by the relative velocity W. The triangle can then be checked by following the 
relative velocity vector and the blade speed to ensure that the same point is arrived at if V alone was 
followed. 


The blade therefore experiences a velocity of magnitude W and angle 8 which will produce a 
force on the blade. 


Recall that from basic mechanics that a force in two dimensions can be resolved into two per- 
pendicular components of any orientation. Two particularly convenient directions are found to be 
perpendicular and parallel to the incoming flow. The force perpendicular to the incoming flow is 
known as the lift force L and the force parallel to the incoming flow is known as the drag force D. 
The principle reason that these are useful directions is that a large body of data on aerofoil perfor- 
mance is available in this form. 


For the wind turbine to produce a useful output a force in the tangential direction must be pro- 
duced, so the lift and drag forces must be resolved into the tangential direction to give a tangential 
component of the force Fg. It is also possible to determine the axial force Fy on the blades which is 
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Figure 3.2: Wind Turbine Blade and Velocity Triangle 


Figure 3.3: Forces on a Wind Turbine Blade 


important for determining the loads on the wind turbine tower and wind turbine bearings. 


All of this is shown in Figure 3.3 which shows the lift and drag forces on the blade. Note that 
these forces are perpendicular and parallel to the incoming flow W and not the axial chord line of 
the aerofoil. Some trigonometry shows the relationship between lift and drag and tangential and axial 
force to be: 
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F; = Lsin 8 + D cos 8 (3.2) 


If we consider Fg and Fy to be forces per unit span. We can obtain a torque per unit span from 
the force times radius. 


T = Fy xr 


The power per unit span is given by the torque multiplied by the rotational speed, this can then be 
integrated from the hub to the tip to give the total power on the blade. 


Ti 
P= Forwdr 


Th 


Where rp, and r; are the hub and tip radius respectively. Usually these conditions are evaluated 
at a number of points along the span of the blade and the integral is obtained numerically by using 
the trapezium rule. The tangential force depends on the lift and drag of an aerofoil which is usually a 
very complex function so no analytical solution is possible. 


In order to determine the performance of the wind turbine some method of determining the lift 
and drag forces is required. This largely comes from test data obtained in wind tunnels, in the form 
a lift and drag plot against incidence. To use this a short digression into aerofoil performance is 
required. 


3.1 Aerofoil Operation and Testing 


In wind tunnel testing an aerofoil is placed in wind tunnel and the incidence of the aerofoil is changed, 
usually by rotating the aerofoil. A lift force perpendicular to the incoming flow and a drag force 
parallel to the incoming flow are measured. Figure 3.4 shows such an aerofoil at two incidences one 
of which is zero or aligned with the incoming flow. Note how the lift and drag remain in the same 
direction with changing incidence and that the aerofoil chord c also does not change with incidence. 
Since it is the flow relative to the aerofoil that produces the lift the incoming velocity is the relative 
velocity W and not the absolute velocity V. The incidence i is defined as the angle made between 
the axial chord of the aerofoil and the incoming flow. 


The lift and drag are usually expressed in turns of a non-dimensional coefficient so that they can 
be scaled for size, fluid density and incoming fluid velocity. These are given by: 


L 
CL = —— 3.3 
L P (3.3) 


where L is the lift force per unit length of aerofoil. The drag coefficient is given in a very similar 
form: 


D 


Cp = ———_— 3.4 
D 1 pW2c ( ) 
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Figure 3.4: Aerofoil at Two Incidences 


where D is the drag force per unit length of the aerofoil. For wind turbine analysis and design val- 
ues of Cz and Cp as a function of incidence 7 are required, these are found in reference books such as 
Abbott and von Doenhoff (1959) A simplified example of aerofoil data is shown in Figure 3.5 which 
is fora NACA 0012 aerofoil. NACA was the predecessor of NASA and the four digit designation 
allows the aerofoil geometry to be determined. On-line coordinate generators are available. ! 


Figure 3.5 is a simplified presentation of aerofoil data as there is no indication on the dependence 
on Reynolds number, for serious work primary sources should be consulted. 


An excellent explanation of the basic physics of aerofoil operation is found in Babinsky (2003) 
but for the purposes of this book it is enough to note that there are three areas of interest for this 
aerofoil plot: 


1. Since the NACA 0012 aerofoil is symmetrical when 7 = 0° the flow on both sides of the 
aerofoil follows an identical pattern so the lift is zero. 


2. When the incidence is non-zero but below what is called the stall point the lift coefficient 
increases rapidly and the drag coefficient increases more slowly. In this case streamlines of the 
flow over the aerofoil will still largely follow the geometry of the aerofoil. 


3. When the incidence reaches a certain level the aerofoil stalls. The streamlines over the aerofoil 
no longer follow the geometry of the aerofoil, the boundary layer has separated and there is a 
substantial reduction in lift along with a substantial increase in drag. 


A common design choice is to place the design point at around 80% of the maximum lift to allow 
for some variation in incidence with stalling the aerofoil. 


Example A wind turbine is designed to work at a condition with a wind speed of 10 m/s and an 
air density of 1.22 kg/m. The turbine has blades with a NACA 0012 profile and is rotating at one 


‘http ://www.ppart.de/aerodynamics/profiles/NACA4 .html 
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revolution per second. The blade chord length is 0.5 m. Taking the design point at 85% maximum 
lift condition and ignoring the drag on the aerofoil estimate the power output per unit blade span at a 


radius of 6 m for each blade. 


Solution The maximum lift shown in Figure 3.5 is around 1.3 and 85% of this value is around 1.1. 
Recall that: 


1 
L= CrzpW°e 


So we need to find W the magnitude of the relative velocity. To do this we consult a velocity 
triangle such as the one in Figure 3.2. From this triangle we can see that: 


W=VU274+V2 


U = wr = 1 x 2r x 6 = 37.7 m/s 


V is the wind speed at 10 m/s so: 


W = yv 37.72 + 10? = 39 m/s 
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Figure 3.5: Cz and Cp for a NACA 0012 Aerofoil 
The relative flow angle 8 may also be calculated: 


B = tan”! (+) = tan! (=) = —75.1° 


Note that the angle is negative as the sign convention in this book is that angles are positive in the 
direction of rotation. With the relative velocity calculated the lift force can now be estimated: 


1 1 
L= CrzpW*e =11x 5 x 1.22 x 39? x 0.5 = 510 N/m 


The tangential force per unit span can be calculated: 


Fy = Leos 8 = 510 x cos 75.1° = 131 N/m 


and finally the power output per unit span can be calculated 


P= Fog x wr = 131 x 37.7 = 4939 W/m 


3.2 Wind Turbine Design 


In wind turbine design the turbine designer has three primary variables that he or she can change: 


1. The type of aerofoil to be used. In this book the NACA 0012 aerofoil is used for simplicity 
but there are actually a wind variety of aerofoils available many designed specifically for wind 
turbine use. 


2. The chord length c along the span of the blade. 
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Figure 3.6: Relationship between y and i 


3. The angle that the blade is set to given the symbol y, this is positive in the direction of rotation 
as for the flow angles a and 8. This is also known as the blade twist. 


Cz and Cp are obtained with a stationary aerofoil in a wind tunnel and the incidence 2 is defined 
for convenience in the wind tunnel testing environment. In order to analyse wind turbines some way 
of relating the incidence 2 and the blade twist angle y is required and this is obtained in the velocity 
triangle shown in Figure 3.6. 


Note that the lift force L is perpendicular to the incoming flow W and so is not at right angles to 
the chord line of the aerofoil. From Figure 3.6 it is apparent that: 


(=p =y 


This arrangement often confuses students as the sign conventions for 2 on the one hand and 8 and 
y on the other hand are in the opposite direction. It is much easier to work in terms of absolute values 
and then determine the orientation from the velocity triangle in that case: 


lil = [yl — 15 | (3.5) 


which is much more “intuitive” when compared to Figure 3.6. 


Example For a wind turbine with the same parameters as the previous example. That is: NACA 
0012 aerofoil, air density of 1.22 kg/m®, chord length of 0.5 m, rotational speed of one revolution 
per second. Calculate the wind speed at which the blade will stall at a 6m radius if the blade angle 
(y) remains constant. 


Solution The first step is to determine y. From the previous example at a wind speed of 10 m/s 
the relative flow angle was given by: 
B= —75.1° 
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The operating point was obtained by examining the lift as 85% of maximum where Cr = 1.1 
from Figure 3.5 the incidence i at this point is around 10°. Using equation 3.5 determine the blade 
angle: 


ly| = |i] + |8| = 10 + 75.1 = 85.1° 


From Figure 3.5 we can see that the aerofoil stalls at an incidence of around 14° which means the 
relative flow angle at the stall condition is given by: 


|| = 85.1 — 14 = 71.1° 


The velocity triangle in this example is the same as that used in the previous example so from 
Figure 3.2 we can see that 8 should be negative so for stall 6 = —71.1°. Using the velocity triangle 
gives a relationship for V: 


= —U —37.7 
pe ( V ) tan  tan(—71.1) amas 


This example illustrates why long turbine blades are twisted, as the radius changes so does the 
blade speed so to maintain the same incidence the blade has to be twisted. 


Example Ifthe blade hub section in the example above at a radius of 1.5m is designed to stall at 
the same wind condition, what would be the local blade angle? 
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Solution So at the hub of the wind turbine the blade angle y must be set to stall at V = 12.9 m/s. 
Since the radius of our analysis has changed so will the blade speed: 


U = wr = 1 x 2r x 1.5 = 9.4 m/s 


From the velocity triangle (Figure 3.2): 
—U —9.4 
=] e [0] 
= — | =t ——— | = —36.08 
2 = tan ( ) an ( 12.9 ) 


Iyl = |G] + listaul —> |7| = 36.08 + 14 = 50.08° 


Using equation 3.5: 


Again using the velocity triangle to get the sign convention correct y = —50.1° 


3.3 Turbine Power Control 


Too high a power output from wind turbines is actually very undesirable. Too high a power output 
can over-stress the blades causing structural failure and the generator in the turbine has only a finite 
capacity to absorb power. Two methods are commonly used to control the power output. 


1. Stall Control. The blade is designed mechanically for maximum lift. At higher wind speeds, 
the blades will stall, causing a reduction in lift and hence tangential force. The drag will 
be increased significantly, which contributes to the axial force, so there is a need to ensure 
the structure can endure the increased axial loading. The drawback of this technique is that 
predicting the onset of stall and the flow around a stalled aerofoil is very difficult. 


2. Pitch Control. The blade is provided with an actuated mechanism to vary the blade angle at 
different wind conditions. When the wind speed is too high, the blade angle is altered to reduce 
incidence. This is shown in Figure 3.7 in Case A on the left the wind turbine blade angle 
y is set so that the incidence 7 is large. In Case B on the right the blade angle has been set 
so that the incidence i has been reduced. The inset in the bottom left of Figure 3.7 shows 
the two blade shapes superimposed to highlight the differences between them. In low wind 
conditions pitching can be used to increase the incidence and hence reduce the power output. 
The drawback here is that a complex actuation mechanism has to be provided. 


3.4 Further Reading 


The method outlined is this Chapter is a simple method of calculating wind turbine performance. The 
simple wind turbine analysis assumes that the wind turbine does not alter the flow going through it 
when momentum changes are considered later it will be clear that this is only a first order approx- 
imation. A more accurate technique is the so called “Blade Element Momentum Method” which 
allows for the variation of flow over the blade and is described in much more detail in Manwell et al. 
(2002) along with much more detail on the general operation of wind turbines. Hansen and Butterfield 
(1993) provide an excellent description of some early development of the aerodynamics of wind tur- 
bine blades. 
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Figure 3.7: Schematic Showing Wind Turbine Pitch Control 


3.5 Problems 


1. A two-blade wind turbine is designed to operate at an atmospheric condition (the air density 
can be taken as 1.22 kg/m) with a wind speed of 22 mph. The turbine blades of 10 m length 
are attached to a nacelle of a radius of 1 m. A preliminary blading design is considered by 
using a NACA 0012 profile with a constant chord length of 1.5 m. The rotational speed of the 
turbine is 30 rev/min. The blading design is taken at a condition corresponding to 80% of the 
maximum lift. 


Calculate the blade angle and the power output per unit blade-span at 20%, 50% and 80% 
spanwise sections for each blade. Note that the blade starts at radius of 1 m so the tip radius 


is 11 m. Estimate the total power output of the wind turbine using the results from the three 
spanwise sections and an approximate integration. 


Answers: —53.8°; —72.4°;—80.8°; 1.201 kW/m; 3.75 kW/m; 7.92 kW/m; 86kW 
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Chapter 4 


Different Turbomachines and Their 
Operation 


Now that the concept of an axial flow turbine has been introduced in some detail the geometries and 
mode of operation of a number of other types of turbomachine will be considered and the cascade 
view and velocity triangles applied to those machines. 


There are three ways of classifying turbomachines, the first is by the type of fluid they work on 
and the second is based on the direction that the flow travels through the machine and the third is 
whether or not they deliver power to or extract power from the working fluid. 


The classification by type of fluid splits machines into two categories. Those that work on com- 
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pressible fluids (most often air or steam ) such as aircraft engines, stationary gas turbines, steam 
power plants and high speed fans. The complication with these machines is that if the velocity of the 
fluid goes above the local speed of sound in the fluid shock waves may result. The second category 
is those that work in incompressible fluids such as water or oil such as pumps or hydraulic turbines. 
Shock waves are not an issue in liquids but the complication here is a phenomenon known as cav- 
itation which occurs if the pressure of the fluid falls below a certain pressure. This is discussed in 
Chapter 10. Machines such as wind turbines and low speed fans fall outside this classification system 
as although the fluid is compressible in practise the changes in density are minimal so they can belong 
to either category. 


The classification of machine by flow direction is more complex, essentially if the flow direction 
is roughly along the direction of the axis of the machine the device is called an “axial” machine 
but if the flow direction turns and has a component substantially in the radial direction the machine is 
know variously as a “radial” or “centrifugal” machine. The meaning of this classification will become 
clearer as the different machine types encountered in practise are described. 


The final classification is that some turbomachines absorb power from the fluid (turbines) and 
some deliver energy to the fluid (compressors,fans or pumps). 


4.1 Axial Flow Machines 


A simple axial flow machine was used to introduce the idea of a turbomachine in Chapter | essentially 
an axial flow machine is one in which the fluid remains parallel to the axis of rotation as it passes 
through the machine. There are many examples of this type most aircraft engines use axial flow 
devices as do all large power station equipment and wind turbines are the most visually striking 
example. 


The complication with axial flow machines compared to the simple one see in Figure 1.2 is that 
more than one stage may be present. This idea allows energy to be extracted in “small bites” from the 
machine allowing very large pressure ratios to be used in a single device - indeed it was this concept 
that is the key to the superiority of steam turbines compared to reciprocating engines. A single stage 
is sufficient when the fluid does not have a great deal of energy to be extracted for a turbine or when 
the requirement is for a modest pressure rise such as in most fans. 


4.2 Radial and Centrifugal Flow Machines 


To introduce the idea of a radial flow machine we will discuss a specific instance, that of a radial 
flow pump. A general arrangement drawing of the pump is shown in Figure 4.1, on the left the cross 
section of the pump is shown and on the right a plan view. Two components are shown in each view 
the outer casing also called the “volute” or “scroll casing” the aim of which is to ensure an even 
distribution of flow out of the machine, the area is therefore gradually increased up to the discharge 
pipe. The second component is the rotor also called the “impeller”. 


The operation of the machine is such that flow approaches the device through the inlet pipe in the 
axial direction, the fluid is then turned through 90° to the radial direction where it enters the impeller. 
The rotor then increases the angular momentum of the fluid and it exits in the radial direction into 
the volute. The pump shown in Figure 4.1 has no stator following the rotor, the fluid slows and 
experiences a rise in static pressure simply by the action of the increased cross sectional area of the 
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Impeller 


Volute 


Figure 4.1: Radial Pump 


Figure 4.2: 3D View of the Radial Impeller 


volute as the radius increases (An analysis of this can be found in Chapter 10) for large pressure rises 
it is necessary to put in a stator row to control this process. 


A 3D view is shown in Figure 4.2 of the impeller only with a radial and the axial direction marked 
on, this hopefully puts Figure 4.1 in context. The blades of the impeller in Figure 4.2 are what is 
known as radial blades, that is no part of them operate in the axial-tangential plane. The impeller as 
a whole is therefore sometimes called a “radial impeller”. This is best seen in comparison to what 
is often called a “centrifugal impeller” shown in Figure 4.3 here the blades extend from the radial 
direction until they also have an axial component. 


Although radial machines are most often found in pumping applications, radial turbines are 
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Plan View 


3D View 


Figure 4.3: Centrifugal Impeller 


widely used in small gas turbines such as the auxiliary power unit found on most aeroplanes and 
centrifugal designs and are used exclusively for turbochargers on internal combustion engines. 


4.3 Radial Impellers 


The velocity triangle can be applied to radial turbomachines in an identical manner for the axial flow 
machine used earlier. The difference is that the cascade view becomes slightly more abstract and the 
velocity triangles may be drawn in the radial tangential plane as well as the axial tangential plane. 


First consider the radial impeller shown in Figure 4.1 and Figure 4.2. The cascade view for the 
axial flow turbine can be considered as an “unwrapping” of the turbine so that the rotational motion 
of the blades w was transformed into a linear motion with magnitude wrm. We apply an identical 
process to the radial impeller in Figure 4.4 as we did for the axial machine in Figure 1.6 on page 18. 


Figure 4.4 shows three views, the first is a 3D view of the radial impeller, the second shows a 
sketch of the turbine as viewed from upstream with the blades, axial inlet and blade tips highlighted. 
To form the cascade view we can approximate the real turbine with a tip radius R = 0.15m and a 
spacing between the blades of s with a machine with an infinite radius and the same blade spacing 
(or pitch) of s. For the radial impeller the cascade view is formed by simply examining the sketch 
produced at the bottom of Figure 4.4 which shows the inlet and the exit from the blades. 


There are two differences between the cascade view for a radial machine and the cascade view 
for the axial machines: 


1. As we shall see later in Chapter 5 there are only two radii that we need to examine to determine 
the work input that the pump applies to the fluid: the radius at where the blades start and the 
radius at the blade tip. Unlike for axial machines the cascade view need only be drawn once no 
matter how large the radius of the machine. 
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Figure 4.4: The Cascade View for a Radial Impeller 
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Figure 4.5: Velocity Triangles for a Radial Impeller 
2. The cascade view for radial machines is formed in the radial tangential plane and not the axial 


tangential plane. 


Having developed the cascade view the four step method for drawing velocity triangles (Section 
2.3) is applied to the radial impeller at inlet and exit, this is shown in Figure 4.5 which shows the 


sketch of the cascade view and the two velocity triangles. 
The cascade view is the same as Figure 4.4 simply turned through 90°. The inlet to the rotor is 

station 1 and the exit from the rotor is station 2 there is no stator in this case so no station 3 is drawn. 

A station 3 does exist however at the volute casing and Chapter 10 will show how to calculate the 
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conditions using conservation of angular momentum - for now we stick to the velocity triangles. For 
radial machines the blade velocity U is different at station 1 and station 2. The blade velocity is given 
by U = wr and will get larger with increasing radius at the same revolutions per minute. 


First consider the velocity triangles at station 1. The first step of the four step rule (Section 2.3) 
is to draw the flow that you know. At station 1 the flow that we know is the absolute inlet velocity at 
inlet to the impeller. There is no mechanism for the fluid to have been imparted with any tangential 
velocity (the inlet flow in the pipe is assumed to have no rotation in it) so the flow is purely radial. So 
a purely radial vector V, be drawn. The second step of the four step rule is to draw the blade speed, 
the blades move in the tangential direction which is perpendicular to the radial direction so we can 
draw a vector Uj. Step three is to close the triangle with the remaining vector w. Finally check the 
correct orientation of the triangle by following the line drawn by W then U this should produce the 
same movement as following the vector V. In Figure 4.5 the velocity triangle is correct so there is no 
need to adjust. 


The velocity triangles at station 2 are drawn with exactly the same method. Step 1: the flow that is 
know this time is the relative velocity, in this book we assume that the flow always follows the blade 
geometry directly so once the blade has been sketched the relative velocity vector Wh at exit from the 
impeller may be drawn. Step 2: The blade speed is always purely tangential but note that the blade 
speed at 2 is greater than at station 1. Step 3: This time close the triangle with the absolute velocity 
vector V3. Finally in Step 4 check that the key rule V=U+W applies by following vectors. 


An alternative approach to step three and four which is essentially a trial and error approach to 
getting the velocity triangle correct is to work out which way around the triangle should go using 
vector addition and get it right first time. The author’s experience is that many students find velocity 
triangles very confusing at first and a slower but fail-safe method is more appropriate when first 
introduced to the concept. 
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Figure 4.6: Common errors in Velocity Triangles 


Figure 4.6 shows correct and incorrect velocity triangles for the radial impeller exit. On the 
left that are two correct triangles the difference being in whether the blade speed is above or below 
the other vectors. Note how the angles and lengths in both triangles are the same. The right hand 
side of Figure 4.6 shows a number of velocity triangles that are incorrect. The first has the relative 
and absolute velocity interchanged, the second two have the vectors for relative velocity and blade 
speed c connected i in an incorrect manner. The incorrect velocity triangles do not obey the key rule: 
V =U + W a tool which gives a systematic method for checking velocity triangles. 


4.4 Centrifugal Impellers 


We now apply the same techniques for deriving velocity triangles to the most complex case commonly 
encountered in turbomachinery - that of a turbomachine where the flow is turned from axial to radial 
and the blades operate in both the axial tangential and the radial tangential plane. The velocity triangle 
at inlet and exit are in different planes. The construction of the cascade view is illustrated in Figure 
4.7 as before. 


In the top left of Figure 4.7 the 3D view of the centrifugal impeller is shown. The front on 
“sketch” is show on the top right, hopefully the relationship between the two is clear as since a 3D 
sketch is hard to do a front view of the 3D model is used. To draw the cascade view two views in the 
radial-tangential plane are shown at the bottom of the figure. 


The bottom left of Figure 4.7 shows a sketch of the blades when they are unwrapped to an in- 
finite radius. The blade tip view is straightforward as it is the same as the blade tip view for the 
radial machine, the complication is that since the geometry is complex the section at the tip only is 
sketched. ! 


The blade hub view (bottom right of Figure 4.7) is actually very similar to the one found in Figure 


!Chapter 5 will show to determine the performance of a turbomachine we only need information at the inlet and exit of 
a set of blades. 
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Figure 4.7: Constructing the Cascade View for a Centrifugal Impeller 


1.6 except the sketch only covers the front part of the blades, the parts that are directly facing us in 
the top right view in Figure 4.7. It is apparent that since the blades extend outwards in the radial 
direction that there will be some variation in conditions in the radial direction at inlet to the blades. 
The blade speed wr will increase with increasing radius. Normally centrifugal impellers are small 
enough to use a mid-height radius as representative of the flow at inlet but for detailed design work 
or for machines with very large blades a number of different radii should be examined at inlet. The 
cascade view is formed (as for an axial flow turbine) by looking down from the casing to the hub so 
you get a plan view of the blades at the inlet to the impeller. 


The two different cascade views are produced in Figure 4.8 along the corresponding velocity 
triangles at inlet and exit. The triangles are built up using the four step rule, the process for which is 
the same as for radial machines. 


For station 1 the cascade view is sketched, it is important to remember that the cascade view is 
only for the entry to the machine, the complex geometry that extends beyond the start of the rotor 
row blading is ignored. The cascade view at inlet is sketched in the axial-tangential plane. At station 
1 the known flow is the absolute velocity which with the absence of any device to increase angular 
momentum will have zero tangential velocity. The blade speed is given by wr, where rı is mean 
radius at station 1 as discussed earlier. With the blade speed and absolute the relative velocity can be 
used to close the triangle and the velocity triangle is complete. 


For station 2 the view is the axial-radial plane and the known velocity is the relative velocity 
leaving the blades. The blades in the centrifugal impeller shown are radial at exit or to put it another 
way the blade direction at exit is parallel to the radial axis. Since the blades are radial the relative 
velocity at exit is also radial. The blade speed is given by wr2 which since rg > rı will be greater 
than the blade speed at inlet. The relative velocity and blade speed allow to complete the velocity 
triangle. 


Figure 4.8 represents the most complex velocity triangle and cascade view encountered in prac- 
tise so once the reader has mastered this technique you are well equipped to analyse almost any 
turbomachine. 
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Figure 4.8: Velocity Triangles for a Centrifugal Impeller 


Note that the analysis of radial flow machines is not restricted to pumps, the technique is equally 
applicable to hydraulic turbines, turbochargers and so on. 


4.5 Hydraulic Turbines 


Hydraulic turbines are devices for extracting energy from a reservoir of fluid usually obtained from 
rainfall above sea level. Water is collected in a reservoir at a height then is carried by a series of pipes 
to a turbine which discharges into a river (Figure 4.9). Hydraulic turbines are universally used to 
generate electricity and there are broadly three types of hydraulic machine in use: The Pelton wheel, 
Francis turbine and the Kaplan turbine, each of which were named after their inventors are are in 
common use around the world today. 


Using the relationship Ap = Ahpg it is possible to express any pressure in terms of height of a 
particular fluid. This is an approach commonly used in hydraulic turbines as pressures directly relate 
to the height of the reservoir or the height above river level that the turbine is situated at. 


Hydraulic turbines are usually classified in terms of total head. For any position in a fluid system 
the total head is given by the following equation: 


H=—+_—+2 (4.1) 


If we examine the total head at the dam of a reservoir (Station 0 in Figure 4.9) where the height 
datum is the river level at exit from the turbine. The surface area is very large so V ~ 0 and since 
the pressure is at atmospheric pressure p = 0. The total head is therefore given entirely by the height 
term z. At entry to the turbine however the height term z ~ 0 but the pressure and velocity terms will 
be much larger - essentially gravitational potential energy has been traded for pressure and velocity. 
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Figure 4.9: Schematic of Hydro-Electric Scheme 
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Figure 4.10: The Four Major Types of Hydraulic Turbine 


Assuming that there are no losses the maximum velocity that would result from a nozzle placed 
next to the turbine open to the atmosphere can easily be found: 


V2 V2 
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pg 2g pg 2g 


But z1 ~ 0 and pp = pı = 0 and vo © 0 therefore: 


V2 
zo = >> Vi = y 2920 


2g 
The four major types of turbines are shown schematically in Figure 4.10. 


The Pelton wheel is the simplest form of water turbine and is used for high head installations. The 
Pelton wheel consists of one or more nozzles which produce a high velocity jet of water these then 
impact a series of buckets which divide the flow in two and the flow exits to the side of the machine. 
The original patent from Pelton in 1880 is shown in Figure 4.11 which includes details of his bucket 
design. The analysis model of the blade is shown on the right hand side of Figure 4.11 essentially 
the blade experiences a relative velocity W1 onto the blades which is then turned through some angle 
and leaves with velocity KW, where k is an empirical coefficient for the reduction in velocity due 
to friction with the bucket etc. Analysis of the Pelton wheel is most convenient when the angle O is 
used in the analysis. The flow is turned through 180° — © and so the Pelton wheel is a high turning 
device as typical values of © are ten to twenty degrees. 


The Francis turbine is used for moderate head installations and in many ways can be though of as 
a centrifugal impeller working in reverse. A cross-section of the device is in Figure 4.10 flow enters 
through a volute or scroll casing which is designed to evenly distribute the flow around the periphery 
of the inlet guide vanes. The inlet guide vanes increase the angular momentum of the fluid which is 
reduced again in the turbine rotor that turns the flow from the radial to the axial direction. The flow 
then exits into a draft tube and from there to a river. 
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Figure 4.12 shows a 3D CAD drawing of a Francis turbine, the real machine is quite complex so 
an external view and a view with the top half of the turbine removed are shown. The Francis turbine 
shown has two sets of guide vanes the first at a high radius are fixed in place, the second allow control 
of the machine by being pivoted on an axis. Usually these guide vanes can turn to 90° so can actually 
stop all flow through the machine to allow it to started and stopped. The control mechanism for these 
vanes is clearly seen in the external view. The relationship between the three dimensional object and 
the sketch in Figure 4.10 should be clear. 


The Kaplan turbine is used for low head applications, the operation is very similar to that of a 
Francis turbine but the rotor design is different. The blades do not extend to the radial direction, so 
the flow is turned before it enters the rotor blade row. This can be seen in Figure 4.13. The Kaplan 
turbine impeller is a lot like that of a propeller. The bulb turbine is very similar to an axial flow 
turbine, except that the medium is water rather than a gas. Again both turbine types are shown in 
Figure 4.10. 


4.6 Common Design Choices 


One common design choice is to have constant axial velocity: Viz = Vox = V3x 


Hence the blade height for turbine must increase as density falls, as seen earlier. There is no 
fundamental physical law that produces this constraint it is a design choice, although it is often used. 


For multi-stage machines, successive stages are often designed to have identical flow angles and 
velocities at corresponding positions in each stage: 


V= V3 ay = a3 and Vig = V39 
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Figure 4.11: Pelton’s Patent Application and Analysis Model 
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Figure 4.12: Three Dimensional Views of a Francis Turbine 
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Figure 4.13: Three Dimensional Views of a Kaplan Turbine 


This is know as a Repeating Stage where the velocity at inlet is equal to the velocity at exit. 


Often we design for axial leaving velocity from a turbomachine. Since for a given V, the exit 
velocity V3 is given by V3 = ,/V2+ Va the minimum exit kinetic energy at exit occurs when 


Voz = 0 or the flow is entirely in the axial direction. Clearly this condition occurs at a3 = 0 so an 
axial leaving condition specifies the shape of the exit velocity triangle for the stator or rotor row. 


These three commonly used design conditions are design choices, although there are many suc- 
cessful turbomachines in use that do use these design choices there are some that do not. 


4.7 The Turbomachine and System 


The turbomachine always operates as part of a system. Pumps are required to deliver fluid at a higher 
pressure for water delivery for example, steam turbines receive steam from a boiler and deliver it to a 
condenser or some industrial process. The torque from a turbine is either used to generate electricity 
or provide mechanical power for gearboxes. The mechanical work for a compressor or a fan comes 
from either electrical power (or in the case of a gas turbine engine is taken directly from the turbine). 
So both the fluid movement and the work changes from the machine fit into a system. 


The apparent exception to this rule are devices such as wind turbines which although they supply 
electricity to some sort of load have varying and unpredictable inlet conditions. Even here a great 
deal of work is done to ensure that turbines are sited appropriately so that they extract the maximum 
amount of energy at a given location and have a minimal visual impact. 


The inputs and outputs from the turbomachine are not arbitrary and are usually fixed by external 
parameters: 


e For steam turbines the inlet conditions are fixed by the boiler conditions which are dependant 
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on the type of fuel used in the cycle and the overall thermodynamic cycle that is chosen for the 
site. In a steam cycle the exit conditions are fixed by the temperature at which the condenser 
cooling fluid can operate at. 


e For aircraft engines the inlet conditions are determined by the operating envelope of the aero- 
plane and the exit conditions largely determined by the amount of thrust required from the 
engine. 


e For pumps the exit conditions are set by the flow rate required and pressure loss introduced by 
the piping system to which the pump is connected to. In the case of moving the fluid from one 
location to another the pump also has to supply a pressure equivalent to the change in height 
between the two locations as well as any losses in the system. 


e For wind turbines the inlet conditions are set by the location of the wind turbine and the variabil- 
ity of the weather. Of all turbomachines wind turbines have the most variable inlet conditions 
and ensuring that the design performs over a wind range of inlet conditions is key to producing 
a successful turbine. 


e For hydro-electric power systems the available pressure or head is determined by the height 
difference between the reservoir and the river into which the fluid is discharged. (See Figure 
4.9) The flow-rate is determined by the rainfall for the region in which the hydro system is 
located and the capacity of the reservoir. A common strategy is to use the reservoir as an 
energy storage system and run the turbines when electrical demand is high, turn them off during 
periods of low demand and allow the reservoir level to rise. 


How the turbomachine fits into the system is generally beyond the scope of this book but it is 
worth remembering the impact that the system will have on the design of the turbomachine. 


4.8 Problems 


1. Sketch the cascade and meridional views for a centrifugal and a radial impeller. 


2. For the Francis turbine in Figure 4.12 sketch the cascade for inlet and exit. Sketch also the 
meridional view. 


3. For the Kaplan turbine in Figure 4.13 sketch the cascade for inlet and exit. Sketch also the 
meridional view. 
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Chapter 5 


Application of The Equations of Fluid 
Motion 


Engineering study of fluid dynamics (and turbomachinery) is built on three sets of equations the 
conservation of mass (continuity), conservation of momentum and the conservation of energy. 


5.1 Conservation of Mass 


Conservation of mass can be simply expressed as “what goes in must come out”, in steady state 
operation the mass flow of the fluid entering the machine must be matched by mass flow exiting the 
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machine. 
In equation form this is expressed as: 


th = pAV (5.1) 


Where rm is the mass flow, p is the fluid density and A is the cross sectional area and V is the 
velocity. For the equation to be valid the velocities and cross sectional areas must be perpendicular 
to each other. This will become clearer with the application of the equation to various machines. 


The simplest case to examine is that of an axial flow turbine. Consider the meridional view of 
the machine shown on the right hand side of Figure 2.4. We denote (as before) the conditions at each 
station by a subscript so p2 is the density at station 2 for example. Conservation of mass tells us that: 


MA Vie = p2A2Vor = p3 A3 V3x (5.2) 


Note that the axial component of velocity is used, as this is perpendicular to the cross sectional 
area under consideration. The cross sectional area of an axial flow turbine is given by the area of an 
annulus between the hub radius r, and the tip radius rz: 


A= n(r? — r2) (5.3) 


In turbomachinery applications it is more useful to express the area in terms of the mean radius 
Tm and the blade height b. 


From Figure 2.4 we can see that: 


tm = ne = 2rm =Te+Th (5.4) 


and 
b= Tt —Th (5.5) 


We can rewrite the area (Equation 5.3) and substitute using Equations 5.4 and 5.5: 


A =n(r? — r?) = a(r¢e + rn)(re — ra) = T(2rm)(b) = Armd 


There are therefore two alternative expressions for the cross sectional area in an axial flow turbine: 


A = T(r? — r?) = 2armbd (5.6) 


Conservation of mass provides us with a powerful tool to explain the geometry of many turboma- 
chines. 


Example An industrial turbine operates at an 8.8:1 pressure ratio and a mass flow of 77 kg/s using 
air as the working fluid. The exhaust temperature is at 43° C and the inlet temperature to the machine 
is around 1000° C. The mean blade radius is 0.4 m. The machine is to be designed for a constant 
axial velocity of 200 m/s. Estimate the blade heights at entry and exit of the turbine. 
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Solution This problem requires us to apply the continuity equation at inlet and exit from the tur- 
bine, to do this we need to know the density of the fluid at inlet and exit. This is obtained from basic 
fluid dynamics: 


- = RT (5.7) 
p 


where p is pressure, R is the gas constant and T is the temperature. For air R = 287 J/kgK isa 
good approximation. The temperature has to be expressed in Kelvin. 


At entry to the turbine: 
pı = 8.8 bar = 8.8 x 10° Pa and Tı = 1000 + 273 = 1273 K 
At exit from the turbine: 


p2 = 1.0 bar = 1.0 x 10° Pa and To = 473 + 273 = 730 K 


The corresponding densities are therefore: 


Pp 8.8 x 10° 
PL RT, 287 x 1273 


= 2.41 kg/m? 


p2 _ 1.0x 10° 3 
= zs = 0.49 k 
P2 = RT, 287x 710 am 


Finally we apply the continuity equation to obtain the required blade heights: 


m = pÅVz = p2armbV, =>> b= Ta 
PTT m Vx 


Since mass flow conservation applies to the turbine (what goes in must come out!). This can then be 
applied to the entry and exit of the turbine. 


B 77 
2,41 x 27 x 0.4 x 200 


-_ 77 
0.49 x 27 x 0.4 x 200 


by = 0.06 m 


bg =0.31m 


This is shown approximately to scale in Figure 5.1. The turbine is split into three stages each with 
a stator and rotor. The design is for a constant mean radius although it is equally possible to design 
for a constant inner or outer diameter. The key point is that Equation 5.1 can explain why the blade 
height changes through a typical gas turbine. This is based on the assumption that the axial velocity 
remains constant, a common design choice that is discussed later. 


5.1.1 Application to Radial Machines 


The application to radial machines is straightforward but requires a good understanding of the geom- 
etry of these machines. Figure 5.2 shows the meridional view of a radial and a centrifugal impeller. 
For the radial turbine (shown on the left) the application of continuity is straightforward. The flow 
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Axis of rotation 


Figure 5.1: Meridional View of a Gas Turbine 


perpendicular to the cross sectional area of interest is the radial velocity and the cross sectional area 
is given by the area of a cylinder face with a depth given by the blade height b and a circumference 
given by 27r, that is A = 27rb so the continuity equation becomes: 


PAVir = p2A2Vor => p12nr1b1Vır = p22nr2b2 Vər 


So for machines operating at constant density and constant radial velocity then bz < bı. The 
blade height falls as the radius increases. 


The application to centrifugal machines (right hand side of Figure 5.2) requires slightly more care. 
At exit the continuity equation is the same but at inlet different velocities and areas are needed. At 
inlet the velocity is the axial velocity and the perpendicular area is the annulus between the hub and 
the tip of the blades. So A = 277,,b where rm is the mean radius at inlet to the machine. Therefore 
the continuity equation becomes: 


my = p1A1 Viz = P1277 mbViz 


and for station 2: 
m2 = p2A2Vir = p22Tnr2bV2ər 


The most useful application of these equations is that given a mass flow and a machine geometry 
the velocities inside the machine can be calculated. 


5.2 Conservation of Momentum 


Newton’s second law of motion can be applied to rotational as well as linear systems in words: “torque 
is equal to the rate of change of angular momentum”. In equation form this is given by: 
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Figure 5.2: Meridional Views of Radial and Centrifugal Machines 
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PE 


(5.8) 


Where T is the torque, L is the angular momentum and t is time. In linear motion force is equal 
to the rate of change of momentum in the system which is the well known form of Newton’s second 
law. 


We now apply Equation 5.8 to the flow through a turbomachinery passage shown in Figure 5.3 
with an entry station 1 and an exit station 2. The passage varies in radius from station 1 to station 2 
and can be any shape (though in the figure the radius change is linear), the passage is bounded by a 
hub and a casing so mass flow is conserved through the turbomachinery passage. 


The velocity through this passage may differ between entry and exit but at either of those stations 
is uniform. The left hand side of Figure 5.3 shows a cross sectional view of the turbomachinery 
passage at entry and at the top a representative velocity vector decomposed into radial and tangential 
components is shown. This vector would have the same value of radial and tangential components no 
matter where on the circumference of the passage at entry we chose to measure it. This is illustrated 
by drawing the same vector at 3 o’clock, 6 o’clock and 9 o’clock on the cross section at entry to the 
annulus. A similar exercise can be carried out at exit from the turbomachinery passage. 


Angular momentum is given by the moment of momentum or: 


L = mVọr (5.9) 


Where m represents the mass. So at entry to the turbomachinery passage the angular momentum 
is mı Viorı and at exit from the passage is m2V2gr2. So the change in angular momentum between 
entry and exit is: 


mə2Və9r2 — Mı Vigri 
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and the time rate of change is given by: 


Mm2Vor2 — Mı Vigri 
t 


= m2Vagr2 — Ma Vigri = m(Vaer2 — Vier1) 


Since “what goes in must come out” mı = mg. The time rate of angular momentum is equal to 
the torque we have: 
T = m(Vaer2 — Vier1) (5.10) 


Recall from basic mechanics that power is torque times rotational speed then: 
P = Tw = mhw(Voere — Vier1) (5.11) 
Where w is the rotational speed of the device. It is very useful in turbomachinery to work in terms 


of quantities per unit of mass going through the turbine, so we divide Equation 5.11 by the mass to 
give one of the most important equations in the book: 


w = w(Voere — Vier1) (5.12) 


This is known as the Euler Turbomachinery Equation. The Euler equation applies to all types 
of turbomachinery, from wind turbines to pumps to gas turbines, whether they be axial flow, radial 
flow or a mixed flow device. It also applies to either energy extraction devices (turbines) or energy 
delivery devices (compressors, fans or pumps). 


There are a number of important points that this equation brings to light: 
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Figure 5.3: A Generic Turbomachinery Flow Passage 


e For a turbine the Euler equation will give a negative value and for a compressor it will give 
a positive value. The energy flows are in different directions so it is reasonable that they be 
of different signs! It is common when dealing with turbines to multiply the right hand side of 
Equation 5.12 by -1 to avoid continually have to write negative signs in front of every calcu- 
lation. Some authors (e.g. Dixon (2005)) refer to different forms of the Euler equation to deal 
with this circumstance but in this book we will stick with Equation 5.12. 


e The power input or output is given entirely by the change in angular momentum. This is 
determined by the turning done by the fluid in the turbomachinery passage which is controlled 
by the blade angles of the turbomachine. In short: work input or output is determined by 
the flow turning - this is one of the fundamental concepts of turbomachinery and is vital to 
understanding the operation of all these important devices. 


e Finally we note that since angular momentum (Equation 5.9) is the product of mass, radius and 
tangential velocity that changes in radius can also influence the power output or input. So radial 
flow machines tend to deliver or absorb much more power per stage than a corresponding axial 
flow device. 


Example A turbine stage with a rotational speed of 3000 rpm is to be designed with an absolute 
inlet angle of 60° and an absolute exit angle of —60° at a mean radius of 0.4 m. The machine is to be 
designed for a constant axial velocity of 450 m/s. Estimate the specific work from this stage. 


Solution This is a straightforward application of the Euler work equation (Equation 5.12), but 
since the mean radius is constant r4 = r2 and we can simplify somewhat. 


w = w(Voer2e — Vier) = wr(Vo9 — Vio) = U(V29 — Vie) 


Here the blade speed is obtained easily: 


2 
U = 3000 x a x 0.4 = 125.6 m/s 
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Lift 
ift 


Figure 5.4: Isolated Aerofoil compared to a Cascade 


The tangential velocities are obtained from the velocity triangle. This example is straight forward 
the absolute inlet and exit angles are given. From Chapter 2: 


Vo = Vz tana 


So at inlet to the stage: 
Vig = 450 tan 60 = 779 m/s 


and at exit from the stage: 
V29 = 450 tan —60 = —779 m/s 


The specific work from this stage is given by: 


w = 125.6 x (—779 — 779) = —195.8 kJ/kg 


5.2.1 The Difference Between a Single Aerofoil and a Cascade of Blades 


This is a good point to explain the differences between a single aerofoil and the operation of a linear 
cascade. The operation of a single aerofoil is much closer to most common experience - nearly every- 
one has been on an aeroplane for example. In this book we do not discuss the detailed mechanisms 
of how an aerofoil works there are other excellent explanations of this phenomenon such as Babinsky 
(2003). Suffice to say that when an isolated aerofoil (Left hand side of Figure 5.4) meets an oncom- 
ing stream of fluid a lift force is produced. A small drag force is also produced. The aerofoil exerts 
and equal and opposite force onto the stream of fluid but since the mass of the incoming stream is 
effectively infinite the net deflection of the flow is zero. 


In the case of the cascade shown on the right hand side of Figure 5.4 the aerofoil operates on a 
finite mass of air as the upper and lower blades of the cascade limit the area that a single blade can 
operate on. With a finite mass the reaction to the lift force produces a net deflection on the flow. 
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This actually provides an alternative derivation of the Euler Turbomachinery Equation. Consider 
the control volume drawn around the blade on the right hand side of Figure 5.4. The force exerted 
on the fluid is given by the rate change of momentum of the fluid - which is Newton’s second law 
of motion. The momentum at stage | is mV and the momentum at stage 2 is moəV2, since we are 
normally interested in the force in the tangential direction and the velocities are normally steady the 
rate of change of momentum and force in the tangential direction is given by: 


Fo = ma Va z maiki m(Vo9 — Vio) 


Torque is then given by the tangential force times the radius and power is the torque times the 
rotational speed: 
T = Form => P= Formw = rmwm(Vr — Vie) 


Again we divide through by the mass flow to give the power per unit of mass flowing through the 
turbomachine: 


P 
w = — = Tmw(Vo6 — Vie) (5.13) 
m 


Equation 5.13 is a less general form of the equation obtained earlier using angular momentum 
(Equation 5.12) as it relies on the assumption of constant radius between entry and exit of the control 
volume - something that is not true of all turbomachine geometries. 


5.3 Conservation of Energy and Rothalpy 


The conservation of energy can be applied to a turbomachinery blade row as well, this leads to a 
new concept called rothalpy. For simple design a very large amount of information can be obtained 
without recourse to Rothalpy but the concept simplifies more advanced work. 


The steady flow energy equation is applied to the flow between station 1 and station 2 of a turbine 
or a compressor. The blade may or may not be in a rotor row. The steady flow energy equation applied 
to the turbine shown in Figure 2.5 is given by: 


"UvA tlasa) (5.14) 


ġ +w =m | (h2 m)+5 


There are two simplifications that can be immediately applied to this equation. Firstly any height 
changes through the machine will be negligible compared to the changes in velocity and enthalpy so 
the z2 — zı term will be zero. 


Secondly the heat transfer to and from the fluid in the machine will also be zero, ġ = 0. This 
second assumption requires some explanation as we also apply this to turbomachines where the tem- 
perature of the fluid is around 1500 K. The reasoning is that the velocity of the fluid is high enough 
that the amount of time a given particle of fluid spends inside the machine is very small so there is 
no opportunity for heat transfer for take place. In reality of course the outside of a high temperature 
turbomachine will get very hot but the heat transfer through the machine casing will be very small 
compared to the work transfers taking place inside the device. 
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Equation 5.14 therefore becomes: 


tb = tin | (ha — h) + F02 va) =r | (m+) ¢ | 7) 


If a moving fluid is brought to rest in an isentropic (or fully reversible) manner different values 
of temperature and pressure will be obtained than if measurements were taken with the fluid still 
moving. The values of temperature and pressure thus obtained are known as the stagnation conditions 
and those obtained whilst the fluid is moving are known as the static conditions. The same concept 
of static and stagnation conditions can be applied to the quantity enthalpy: 


y2 
ho = h + T (5.15) 
where ho is known as the stagnation enthalpy and A is called the static enthalpy. 


Equation 5.14 can be further simplified to: 
w =m [hoz — hor] (5.16) 
Equation 5.16 is important as it relates the work output or input from the turbomachine to the 


change in fluid properties between the inlet and the exit of the machine. The method for determining 
ho varies depending on the type of fluid. 
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5.3.1 Rothalpy 


Note that w/7n is the specific work w so the equation is further simplified to: 
w = hoz — hor 


But the Euler Turbomachinery Equation (Equation 5.12) also contains an expression for w, so 
Equation 5.14 becomes: 


w = UgVo9 — Ui Vie = hoz — hoy => hoi — UiVie = hoz — U2V26 


Define a new quantity called rothalpy given by: 


I = ho — UV (5.17) 


From the re-arrangement of the SFEE energy equation, rothalpy is constant from station 1 to 
station 2, I1 = Ig, which leads to the general principle that the conservation of energy applied to 
turbomachines can be expressed as “over a blade row rothalpy is constant”. 


Rothalpy is conserved over a blade row (stator or rotor) and not a stage (stator plus rotor). The 
value of rothalpy at the stator/rotor interface will clearly change as the stagnation enthalpy ho will be 
the same but the blade speed U will be different depending on whether or not one is considering the 
rotor or the stator. 


5.3.2 Rothalpy in Stators and Rotors 


For a stator Rothalpy conservation is straightforward, since U = 0, rothalpy conservation from stator 
entry at station 1 to stator exit at station 2 yields: 


y= hor = hoz 


Or in words - in a stator the stagnation enthalpy is conserved. For a rotor the situation is more 
complex. Consider a generic velocity triangle in Figure 5.5 from which we will derive some general 
relationships between relative and absolute velocity. Since the relationships are general no subscripts 
are used to indicated a particular station, the relationships could apply to any station. Some of these 
were discussed in Chapter 2 when the concept of velocity triangles was introduced but are worth 
repeating here. 


Two relationships are obvious from the triangle: 


Vz = W, and Vg = U + Wọ 


= 


Pythagoras’s rule can be applied to the absolute velocity, V: 


V? = V2 4+ VŽ 
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Figure 5.5: Generic Velocity Triangle 


We can then substitute the two obvious relationships so: 


V? = V? + V? =W2 + (Wọ + U)? = W? + 2WọU + U? 


So we can express the absolute velocity magnitude V in terms of relative velocities and the blade 
speed: 
V? = W? + 2WọU + U? (5.18) 


Now consider a rotor row with entry (point 1) and exit (point 2): 


h =h = ho — UVi9 = hoz — UV26 


expanding out the ho term yields: 


2 2 


V, V. 
hy + -> — UVis = ha + > — UVa 


Consider the left hand side of this equation 


V2 
LHS = hı + a — U Vio 


Substitute for V using Equation 5.18 and for Vg using Vg = U + Wọ: 


wW? U? 
LHS = h, + — -2 
1+ 5 5 


We can simplify this further by defining another new term called the relative stagnation enthalpy 
given by: 
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w2 
horel =h+ D (5.19) 


Apply this to the right hand side of the equation as well as the left rothalpy conservation for a 
rotor row reduces to: 


2 


U2 U. 
horeli — oi = horel? — a (5.20) 


where station 1 is the entry to the rotor row and station 2 is the exit to the rotor row. 


Equation 5.20 brings out two significant points about turbomachinery: 


Firstly for axial machines with little change of radius Uj; = Uy so rothalpy conservation is 
equivalent to conservation of relative stagnation enthalpy. Compare this with the stator case 
where rothalpy conservation is equivalent to conservation of stagnation enthalpy. Constant 
relative stagnation enthalpy implies that no work is done in the relative rotor coordinate system 
- which may be confusing but there is no movement in this coordinate system as the observer is 
attached to the rotor blades - therefore in this coordinate system there is no work. The energy 
inputs or outputs depend on the coordinate system chosen but the pressure and temperature are 
independent of the coordinate system chosen. 


Secondly the U? term can be rewritten: 


and can be thought of as the potential energy of a centrifugal force field due to rotation. With 
the same rotational speed a higher radius object will have a higher tangential velocity. 


5.4 Problems 


1. An industrial gas turbine operates at an 8.8:1 pressure ratio and a mass flow of 77 kg/s. The 
exhaust temperature is at 437°C and the inlet temperature to the machine is around 1000°C. 
The machine is to be designed for a constant axial velocity of 200 m/s. (This data is based on 
the Rolls-Royce Avon, an engine that dates from the 1950s but is still used as low efficiency 
high reliability engine for stationary power such as pipeline pumping.) 


For this gas turbine, draw to scale the meridional view with: 


(a) a constant mean radius of 0.4 m. 


(b) aconstant outer diameter (to minimise cross sectional area) of 1.05 m. 


2. Consider an axial compressor rotor blade row. The inlet velocity is 150 m/s and is in the 
axial direction. The blade speed is 180 m/s and the relative outlet angle is —30o to the axial. 
The axial velocity is constant across the row. Calculate the relative flow angle at inlet and the 
absolute swirl velocity at exit and the absolute flow angle. Calculate also the specific work 
output. Answers: —50.2°; 93.4 m/s, +31.9°, 16.8 kJ/kg 
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Chapter 6 


Efficiency and Reaction 


6.1 Efficiency 


There are over one hundred different ways of defining the efficiency of a turbomachine, though only 
one of them will be discussed in this book. Consider a turbine on a h-s diagram, Figure 6.1. 


Three points are shown on this diagram 1 the entry condition to the machine, 3 the real exit 
condition and 3s the isentropic (or ideal) exit condition. (We use 3 rather than 2 as we will consider 
the stator/rotor interface as 2 later) Since the real exit condition will contain irreversibility’s due to 
fluid friction, heat transfer over a finite temperature difference and so on the entropy in the real case 
increases meaning the line 1 to 3 is a diagonal. The ideal process is a straight line on the chart as the 
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Figure 6.1: Enthalpy-Entropy Diagram for a Turbine 


change in entropy is zero for an reversible process. 


One of the reasons there are hundreds of definitions of efficiency is that there are various ways 
of accounting for the kinetic of the flow at inlet and outlet. Recall that for each location in a fluid 
machine there is a static condition and a corresponding stagnation condition. The stagnation condition 
is the condition that occurs when the fluid is brought reversibly to rest so the entropy is the same 
between the static and stagnation conditions. For the turbine stage on the h-s diagram (Figure 6.1) at 
each condition 1,3 and 3s there is a corresponding stagnation condition 01,03 and 03s which are also 
plotted. 


Note that the difference between static and stagnation conditions is always a vertical line on the 
h-s chart as it is by definition an isentropic process. hg = h + u by definition, so the vertical length 
of the line is uae 

Efficiency can be defined in terms of actual work output and ideal work output. For a turbine, the 


Isentropic Total-to-Total Efficiency is the ratio: 


Actual Work Output wA ho3 — hor 
Ideal Work Output wr ne ho3s — hor 


(6.1) 


The ideal work output is always greater than the actual work output, the real (non-deal) process 
serves to reduce the power output of the device. Note that the work output will be negative which is a 
matter of sign convention, if the work output for turbines is negative the work input for compressors 
is positive - one has to pick one to be negative! Often however designers deal with only turbines and 
make the work output positive to avoid writing minus signs in front of every calculation. 


Figure 6.2 shows the h-s diagram for a compressor or pump. Since compressors deliver energy to 
the fluid the net result is a rise in pressure. For a compressor the ideal work is always less than the 
actual. The isentropic total to total efficiency for a compressor is given by: 
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YN? 


Figure 6.2: Enthalpy-Entropy Diagram for a Compressor 


Ideal Work Input wy = _ hogs — hor 
Actual Work Input wA "hg hos — hoi 


(6.2) 


Example A compressor stage working on air at 1 bar and 25°C has a work input of 17 kJ/kg 
and an isentropic efficiency of 0.9. The velocity at inlet and exit is the same. Assuming that the air 
properties are unchanged over the stage calculate the pressure output of the stage. 


Solution The work input is the specific work input so and since the velocity at inlet and exit are 
the same V3 = V; the difference between static and stagnation enthalpies is the same. 


WA = hg = hy = 17k J/kg 
Since we can approximate air as a perfect gas (Cp = 1.01 kJ/kgk): 


Ah 17 
CpAT = Ah AT = — T; — Ti = — = 16.83 K 
$ er eae Ey! 
The efficiency for a compressor using a perfect gas can be expressed in terms of temperature rises: 
_ fasta _ Ths —Th 
hai day 


From basic thermodynamics we recall that for an isentropic compression: 


=0.9 = Tz,- T) = 0.9(T3 = Tı) = 15.15 K 


n 


a aa S 
Tss _ (2) i 
Tı pı 


1.4 


04 
) = 1.19 bar 


Now y ~ 1.4 for air at 20° so: 


ps3 _ (298 + 15.15 
po 298 


Since the inlet is at atmospheric pressure 1 bar the pressure output from the stage is 1.19 bar. 
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6.1.1 Using Efficiency 


Actual efficiencies can be determined by either measurement or computationally intensive simulations 
of the fluid flow through a turbomachine. For the purposes of this book they will simply be based on 
given values. Isentropic efficiencies of over 90% have been reported for turbines and compressors are 
typically a few percent less efficient. 


The key point is that efficiency provides a link between the work from flow geometry (Euler work 
equation) and the actual enthalpy or pressure change in the fluid. 


For a steam turbine the ideal work is given by: wr = hoi — hoss. Finding the ideal work 
involves determining the enthalpies at inlet and exit. At inlet hı is determined from the h-s 
diagram or steam tables given inlet pressure and temperature. The exit condition hgs is derived 
from an isentropic expansion (vertical line on the h-s chart) and a known exit pressure. 


For a gas turbine the ideal work is obtained using the ideal gas laws: 


ay 
Ti Ea 
wr = hoi i h3s = CpTo1 (1 = 7) = CpTo1 ( = (=) ) 
OL Poi 


e For a pump operating on an incompressible fluid (water or low speed gas flow) then the ideal 
work can also be obtained from an equation: 


A 
wr = hor — hogs = = = gAHo 


where AH is the change in total head across the turbine. 
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Example A steam turbine stage operates with inlet conditions of 120bar and 500°, the exit pressure 
is 10 bar. If the efficiency of the stage is 90% calculate the specific work output. 


Solution The strategy here is to look up hı and hg, from a suitable set of Steam Tables such as 
Rogers and Mayhew (1994) and then use the efficiency to get the actual work out. 


From steam tables or a Mollier Chart: 
hı : 120 bar, 500°C = 3350 kJ/kg 


Sketch a line of constant entropy on the Mollier chart from point 1 to the 10 bar line or interpolate 
using tables to determine hgs: 
h3s = 2730 kJ/kg 


rearranging the definition of turbine efficiency (Equation 6.1) the actual work is: 


wa = nit(hgs — h1) = 0.9 x (2730 — 3350) = —558 kJ/kg 


The answer is negative as our sign convention is that energy transfers out of the system are negative, 
but generally such sign conventions are not important in turbomachinery as one normally knows 
whether the device is supposed to extract or deliver energy! 


6.1.2 Other Efficiency Definitions 


In this book only the total to total efficiency (using only ho) in definitions of efficiency is used but 
using static or stagnation conditions in the definitions of efficiency allows the kinetic energy of the 
leaving fluid to be used in different ways: 


1. Where the exit kinetic energy is useful for example for thrust in a jet engine, or subsequent 
stages in a multi-stage machine. In this case the energy is not part of the maximum work 
available. 


2. Where the exit kinetic energy is not useful for example in the last stage of a steam turbine 
exhausting to the condenser, or a hydraulic turbine where the water flows away done the exit 
pipe. An ideal turbine would make use of that exit kinetic energy so it must be included in the 
definition of the maximum work available. 


These different cases explain why other definitions of efficiency may be used. 


6.2 Reaction 


There are three difficult ideas presented in this book. The cascade view, velocity triangles and the 
concept of reaction. Once these three concepts have been mastered they open the door to a much 
wider understanding of turbomachinery which is why this book concentrates on them. The third 
concept to master is that of Reaction which is defined as the ratio of the static enthalpy drop through 
the rotor divided by the static enthalpy drop through the stage. Or in equation form: 
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_ Ahroror 
AhSTAGE 


Recall from basic thermodynamics that: 


(adh ap 
p 


If the process is isentropic ds = 0 then: 
1 1 
dh = -dp => Apx -Ah 
P P 
Most turbomachines are quite efficient so approximating their operation as isentropic is reason- 


able. This means that enthalpy drops can be approximated as pressure drops and gives the widely 
used definition of Reaction: 


AhROTOR _ ApRoror 
AhsTAGE  APsTAGE (6.3) 


Consider the turbine shown in Figure 2.4 on page 25 and in the context of Reaction a number of 
questions arise: How should we assign the pressure drop of the stage between the stator and the rotor? 
What blade angles do we need to set to obtain the desired pressure drops? What levels of reaction are 
desirable? These are questions that we will answer next and in Chapter 8. 


6.3 Reaction on the h — s Diagram 


It is possible to show the level of reaction of a particular turbomachine on the h-s chart and get a good 
understanding of what is going inside a turbomachine as a result. This can get somewhat complex 
so the chart is built up in stages. Firstly consider Figure 6.3 which shows a real process in a turbine 
expanding from pressure pı to pressure p3. The geometry over which this takes place can be seen in 
meridional and cascade view in Figure 2.4 on page 25. So station 1 is the stator inlet, station 2 is the 
rotor/stator interface and station 3 is the rotor exit. The intermediate pressure ps is shown in Figure 
6.3 along with the static and stagnation conditions at 1 and 3. 


We now consider station 2 on the h-s diagram in more detail. For the stator 1 to 2 ho, = hog 
because no work is done in the stator we can then add a point 02 to a h-s diagram showing this. Since 
the difference between stagnation and static enthalpy is given by % “y we can draw the corresponding 
velocity magnitudes on the h-s chart for Vj,V2 and V3. This done in Figure 6.4. 


Now since the stagnation enthalpy is constant: 


Ve Ve V2 v? 
hor = hoz = h+ =h 2+- 5 = h hm 5 


From Figure 6.4: V2 > Vj so there is a static enthalpy drop in the stator. For a turbine the flow is 
accelerated in the stator. 
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The level of static enthalpy or pressure drop is controlled by the velocity V2 for high values of 
Və the pressure drop is greater. V2 is set by turbomachine geometry so we can control the level of 
reaction in the turbine by altering the blade angles and blade height through the machine. 


Finally we consider a turbine where the mean radius is the same so the blade speed U is the 
same at entry to and exit from the rotor. This restricts the generality of our diagram somewhat but 
this condition is true of a great deal of axial flow turbines and provides some useful insights. If the 
blade speed is the same through the turbine the relative stagnation enthalpy will be constant through 
the rotor blade (See Chapter 5). Two new points are added to the h-s diagram 03R and 02R which 
represent the relative stagnation enthalpy at points 2 and 3 respectively, this is shown in Figure 6.5. 
Points 03R and 02R lie above the corresponding static conditions 2 and 3 and the distance between 
them is given by vi from the definition of relative stagnation enthalpy. Figure 6.5 shows that for an 
axial flow turbine the flow is accelerated through a turbine W3 > Wo. 


Reaction is often a cause of some confusion for students. This is because intuitively one might 
expect the power output from a stage to depend on the pressure drop across it - but as we showed with 
Equation 5.12 the work input or output from a turbomachine depends entirely on the flow turning the 
rotor row. The relationship between pressure drop and velocity change can be seen in the enthalpy- 
entropy diagram in Figure 6.5. 


The effect of the real process irreversibilities can also be seen in Figure 6.5, an ideal process is 
a straight line (As = 0). So if the process was ideal it could be drawn on the h-s diagram using a 
vertical line from 01 to the pressure po, the velocity V2 would then be larger than in the real (As > 0) 
case. 
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Figure 6.3: Basic h-s diagram 


Figure 6.4: h-s diagram with ho 
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Figure 6.5: h-s diagram with ho and horel 


6.4 Problems 


Note that some of these problems use the same data as previous chapter problems so if you have kept 
your answers you may save some time. 


1. An axial compressor has a rotor blade row where the inlet velocity is 150 m/s and is in the 
axial direction. The blade speed is 180 m/s and the relative outlet angle is —30° to the axial. 
The axial velocity is constant across the row. 


These rotor blades are followed by a set of stator blades to form a complete stage. The stators 
turn the flow back to the axial direction and the axial velocity is constant across the whole 
stage. If the total-to-total efficiency is 90%, calculate the total pressure at outlet from the stage. 
Assume air properties with inlet stagnation conditions of 20° and 1.0 bar. Answer: 1.187 bar 


2. Calculate the reaction of the stage in the previous question. Answer: 0.74 


3. The flow at exit from a turbine stator row has a velocity of 100 m/s at an angle (a2) of 70° 
to the axial direction. The rotor row is moving with a velocity of 50 m/s. At exit from the 
rotor row the relative flow angle (83) is —60°. The axial velocity is constant across the row. 
Calculate the power output for a flow rate of 4kg/s. Calculate the total pressure drop across the 
stage if the efficiency is 90 % and the fluid density is constant at 1.2 kg/m3. Calculate also the 
static pressure drop. What is the stage reaction? Answers: 20.6kW, 6.88kPa, 6.93kPa,0.154 
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Chapter 7 


Dimensionless Parameters for 
Turbomachinery 


Dimensionless coefficients are a very useful technique to bring out trends between groups of variables 

in fluid mechanics. Further details can be found in basic fluid mechanics texts Massey (1989) but 

essentially the technique works by constructing physically plausible relationships between variables. 

Consider the relationship between static pressure, dynamic head and stagnation pressure for example: 
pV? 


Le Da 


For this relationship to be physically possible pV? must have the same units as pressure, in simple 
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terms one cannot count two oranges and expect to end up with an apple! 


Pressure can be explained as a force per unit area, force can be expressed as a mass times an 
acceleration so the units of pressure will be given by: 


length 


mass x acceleration mass X Fr mass 
area ~ lengthè —s time? x length 
The units of pV? can similarly be expressed: 
Wi eee mass length? mass 
ensi velocity = = 
z á length? time? time? x length 


Therefore the equation is said to be dimensionally correct. Using this principle one can construct 
a number of dimensionless parameters or dimensionless coefficients. For example reaction is a di- 
mensionless coefficient as the units on the top of the expression cancel the units on the bottom of the 
expression. 


Dimensionless parameters provide a number of advantages: 

1. With an appropriate choice of dimensionless parameters the performance of machines can be 
characterised using only a few key variables. 

2. Given data on one size of machine performance can be predicted at different sizes. 


3. Given data on one set of operating conditions behaviour at different operating conditions can 
be predicted 


4. Enable designers to pick a particular machine shape of maximum efficiency. 


Which dimensionless parameters are important depends very much on the application area several 
of which are discussed in turn. 


7.1 Coefficients for Axial Machines 


For axial flow turbines there are four dimensionless coefficients that have been found to be important. 
The first of these is the flow coefficient: 


pE (7.1) 


which is the axial velocity divided by the mean blade speed. For a given blade speed and height 
this gives some idea of the flow through the machine. 


The second of these is the stage loading coefficient or work coefficient: 
Wr Aho 


Sa (7.2) 


which is a measure of the work done in a stage. The other two important coefficients are the 
efficiency 77 and the reaction R. 
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Case Flow Coefficient Stage Loading Coefficient 
Aircraft Engine Compressor 0.4 to 0.70 0.35 to 0.50 
Aircraft Engine HP Turbine 0.5 to 0.65 1.0 to 2.0 
Aircraft Engine LP Turbine 0.9 to 1.0 1.0 to 2.0 


Table 7.1: Typical Values of Y and ® 


Velocity Triangle at 3 Velocity Triangle at 2 


AV=V V30 


Figure 7.1: Velocity triangles for exit and inlet combined 


In general the designer of turbomachinery blading does not have complete freedom to determine 
the dimensionless coefficients that the design will operate at. For example the mean blade speed is 
often fixed by mechanical considerations, or the rotation frequency of the machine is fixed by the 
frequency of the electrical supply and the number of poles in the electrical generator. The stage 
loading is often fixed by the number of stages and the required enthalpy drop in the machine. 


Table 7.1 shows some typical values of the stage loading and flow coefficients for aircraft engines. 
Note the much lower loading for compressors compared to turbines, the reason for this is that com- 
pressors operate in an adverse pressure gradient. That is compressors have a higher pressure at exit 
than they do at inlet and this makes the boundary layers along the blade surfaces much more likely to 
separate. The aerodynamic design is then much more challenging than for a turbine so lower levels 
of loading are applied. 


For an axial flow turbine the mean blade speed is usually constant across a blade row as there 
are no large radius changes in most axial flow turbines. It is impossible to change the blade speed 
w within a blade row. It is therefore often convenient to plot the exit velocity triangle and the inlet 
velocity triangle on top of one another (Figure 7.1. 


Now the Euler equation (Equation 5.12) gives: 


w = w(Voere — Vier1) = UmAVe 
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So substituting this into the expression for work coefficient: 


and remembering the definition of flow coefficient: 


Vz 
Um 
Two things become apparent: 


e Fora given blade speed Um: Y gives the width of the diagram 


e For constant axial velocity V,: ¢ gives the height of the diagram 


The representation in Figure 7.1 is compact and is often used when designing axial flow ma- 
chines, especially when large numbers of stages have to be put together - in this book the emphasis is 
simplicity not compactness so it is not used very much. 


One of the key points to take away from Figure 7.1 is that given a constant blade speed and axial 
velocity changing the dimensionless coefficients has a direct bearing on the flow angles. So for the 
triangles shown in Figure 7.1 if the stage loading coefficient is reduced for example but the flow 
coefficient and blade speed kept the same, so U and V, are constant then the absolute flow angles a1 
and a3 will reduce. There is a direct link between machine geometry and dimensionless coefficients. 


Finally an important simplification with constant blade speed is that since Vg = Um + Wọ then 
AV = AWọ 


Example Given a turbine blade row with constant axial velocity of 150m /s at 5000rpm on a mean 
radius of 0.7 m and an absolute flow angle at exit from the stator of 70°. The turbine operates with 
axial leaving flow and is a repeating stage. Calculate the flow coefficient, stage loading coefficient 
and reaction. 


Solution This problem is a straightforward application of dimensionless coefficients. First calcu- 
late the mean blade speed: 


2 
Um = wr = 5000 x < x 0.7 = 366.5 m/s 


Vz 150 
i Um 366.5 GAD 
i- w Um(V39 — V29) _ Vse — Væ 
JA UZ Um 


Now the turbine row is axial leaving so V39 = 0 from the generic velocity triangle: 
V29 = Vz tan ag = 150 tan 70 = 412.12 m/s 


so the stage loading coefficient is: 
41212 


366.5 


Y = 1.124 
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Reaction is given by R = Ahroror/Ahsrace. The stage enthalpy drop is given by the work 
output as the turbine has repeating stages so V, = V3. 


AhstT AGE = Aho = wW = Um(AVo) = 366.5(412.12) = 151.04 kJ/kg 
From Figure 6.5 we see that the static enthalpy drop across the rotor is given by: 
W2 = W2 W2 — w2 
Meson = 2 De ce ae 20 
2 2 
Since the turbine axial velocity is constant W2, = W3, and so: 
W39 = Wo 
2 

However since we have constant blade speed AW» = A Wọ so: 


412.12? + 0? 
Ahroror = SIT LO = 84.92 kJ/kg 


Which means that the reaction is given by R = 84.92/151.04 = 0.56 


Ahroror = 


7.2 Coefficients for Wind Turbines 


For wind turbines two coefficients are in widespread use, the first Cp or power coefficient is essen- 
tially an efficiency by a different name: 


Power from Turbine P 
Cp = —— = 7.3 
= Power in Wind 5pV3A oe 
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where P is the power output, V is the velocity of the wind and A is the swept area of the machine. 
The second coefficient is the tip speed ratio, which is the ratio of the wind speed to the tangential 


velocity of the tip of the turbine: 
WP tin 


A= 


(7.4) 


Performance of wind turbines is usually presented as Cp vs A curve, this is often called a “power 
curve”. A power curve for a Nordex N80 turbine 'which produces 2.5MW at full power is shown 
in Figure 7.2. The power curve shows a characteristic of many turbomachinery devices, that there is 
distinct peak in efficiency with a lower device effectiveness away from the design point. 


Example Given the power curve shown in Figure 7.2. Calculate the power output of the device at 
10 rpm in a 10 m/s wind given that the turbine is 80m in diameter and the density of air can be taken 
as 1.15kg/m?. 


Solution This problem reduces to finding the tip speed ratio À looking up the power coefficient on 
the graph and then calculating the power output. The tip speed ratio A gives: 


10x2m ., 80 


WT tip 60 2 
V 4 


So look up the power coefficient that corresponds to 10.47 on Figure 7.2 which is approximately 
Cp = 0.37. The area of the turbine is given by A = mrj,, = T x 40? = 5026 m?, the assumption is 
that the area of the nacelle is negligible. Putting this all together: 


P 


en 
P Tpv3A 


1 1 
— P= Cp3pV°A =0.37 x 5 x 1.15 x 10° x 5026 = 1.069 MW 


7.3 Coefficients for Hydraulic Machines 


Hydraulic machines are those turbomachines that operate with liquid (most often water) as the work- 
ing fluid. These are enormously significant devices which produce up to 20% of the worlds electricity 
in large and small hydroelectric plants. Pumps are vital to the infrastructure of water and fuel supply. 


Dimensional analysis of hydraulic machines yields a different set of coefficients, these are the 
flow coefficient: 
Q 


iÜ 
1 ND? 


(7.5) 


where N is the number of revolutions per second made by the machine and D is the machine 
diameter. We give the coefficients the symbol II. The next dimensionless parameter is the head 


coefficient: 
gH 
Tb = yap? a 


where H is the head produced or absorbed by the hydraulic machine and g is the acceleration due 
to gravity. Finally the power coefficient is given by: 


‘http ://www.nordex-online.com/en/products-services/wind-turbines/n80-25-mw/ 
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Figure 7.3: Collapsing Pump Data onto Non-dimensional Curves 


P 


i= Á 
3 pN3D5 


(7.7) 


The power of dimensional analysis can be seen in Figure 7.3 the left hand plot consists of pres- 
sure rise against flow in the pump at different % of the maximum speed and is taken from a pump 
experiment at the University of Durham. The right hand plot presents the same information in non- 
dimensional terms this time using flow coefficient A compared to pressure rise coefficient ape 
the data is exactly the same but reduces the relationship at different speeds to a single line. This shows 
how plotting parameters in terms of dimensionless quantities is a powerful method for reducing the 
complexity and amount of information an engineer needs to understand. 
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Clearly efficiency is a non-dimensional parameter for hydraulic machines and again it varies 
between pumps and turbines as follows: 


Ideal Work Actual Work 


pump = Actual Work’ "turbine Ideal Work 


7.3.1 Specific Speed for Turbines 
Note that for a turbine: 


Hi — P 
Il pQgH 


For a pump the reciprocal applies. So if for a model turbine we obtain test data and plot the dimen- 
sional parameters against each other the performance of larger machines of the same type can be 
predicted. For a particular machine the operating conditions are expressed by values of N, P and H. 
The rotational speed N is normally fixed by the frequency of electricity supply, P and H are set by 
the flow rate and the height drop of the physical location that the hydro-electric scheme is proposed 
for. It would therefore be useful to have a dimensional group which includes N, P and H but not 
D so that it would be independent of the machine size. This can be obtained by manipulating the 
dimensionless groups for a turbine to obtain a new dimensionless coefficient: 


n 


O mP NPI? 
“Mie peA 


IM; 


The point of machine design is to obtain maximum efficiency and there is generally only one set 
of values of IT; and Iz and IIs for which this occurs. For a given shape of machine we are therefore 
interested in a unique value of I4. We call this dimensionless specific speed: 


1/2 
= I! = NPE (7.8) 
a cers pram ; 


Ks, has a unique value for maximum efficiency for a given shape of machine. It is a shape 
parameter independent of the size, D. By calculating the value of K, from the design specifications, 
we can determine the most efficient machine type. 


Normally hydraulic turbines use water and remain on the surface of the earth so p and g are 
dropped to give a dimensional specific speed: 


NpP!/2 
s= -H5 (7.9) 


For N in rev/min, P in kW and H in m, Ng = 1042K,. 


The term specific speed arises from the idea of a given type of machine producing unit power 
(1kW) at unit head (1m), N, is the speed it would run at in rev/min. It is much better to think of it as 
a shape parameter. 
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For commercial machines there are size limitations the upper limit for N is determined by material 
strength (centrifugal forces) and machines with a small power output P of a given type become less 
efficient - giving a range of applicability for different types of hydraulic machine. 


7.3.2 Specific Speed for Pumps 


For a pump the requirement is more likely to be flow Q, rather than power so a different dimensionless 
group is picked. 
i 7 Q13 n2/3 


I; = = 
nm? (gH)? 
raise to the power 2/3 to get: 
NQ! 
Ky = (gHy3/4 (7. 10) 


which is the dimensionless specific speed for a pump, which again has a unique value for maxi- 
mum efficiency. As for turbines there is also a dimensional specific speed: 


NQ! 


Ns = H 


(7.11) 


For N in rev/min, Q in m/s and H in m, N; = 333K y 


Note that the specific speed for pumps and turbines are related by a simple ratio - they are funda- 
mentally the same concept just applied to different machines. 


7.3.3 Using Specific Speeds 


Pump and turbine manufacturers will provide charts and diagrams that relate pump or turbine ge- 
ometry to specific speed. These charts are based on the manufacturer’s experience with particular 
machines as the choice of machine type depends on the maximum efficiency which can only be de- 
termined by testing. Note that different companies or organisations have widely different choices for 
the units they apply to their charts - be sure to use the same units as the person who made the chart! . 


Figure 7.4 shows a plot of specific speed against head from a large number of hydroelectric 
power plants produced by Alstom Power from 1965-2006. The units of specific speed are actually 
kg'/2m/> s—5/2 or in other words not very meaningful so they are not recorded on the graph. The 
graph shows clearly the expected trend of specific speed vs head. This can also be obtained from 
Equation 7.11 where Ns x H~'° or approximately 1/H in order. The point of the graph is to show 
how the different turbine types are grouped together.” 


From this illustration of industry practise one can draw up some guidelines for turbine type based 
on specific speed, these broad guidelines are shown in Table 7.2. 


>The outlier on the graph with a head of 10m and a specific speed of around 450 is actually a tidal stream turbine at La 
Rance in France so operates with the flow going in either direction so the design is unusual. 
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Case Specific Speed 

Pelton Wheel 50 and below 
Francis Turbine 100 to 300 

Kaplan Turbine 400 and above 


Table 7.2: Values of N, for Different Turbine Types. 


Example A turbine is to be designed for a site with 400 m of head an an expected power of 1 MW, 
the turbine will feed electricity into a 50 Hz electrical grid. Using specific speed estimate which sort 
of turbine design should be investigated further. 


Solution The strategy here is to calculate the specific speed. The head and power are given so we 
need to find possible values or rotational speed. Electrical generators typically run at fixed multiples 
of turbine speed the fastest being a two pole generator which runs at the frequency of the grid, there 
is no limit to the number of pole pairs that can be used in the generator but up to ten pole pairs will 
be examined. 


With a single pole pair, the rotational speed will be N = 50 x 60 = 3000 rev/min: 


= NPY? 3000 x 10001/? 


8 = Spe a = 33.08 


Note that power is input in kW, rotational speed in rev/min and head in m. With ten pole pairs, 
the machine rotational speed will be N = 300 rev/min and so: 
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Alstom Major References 1965-2006 


600 E Pelton 


500 ® Francis 
Y Kaplan 


400 Ww å Bulb 


Specific Speed / [units] 
> 
<j 


a Oi 
10 100 1000 10000 
Head / Im] 


Figure 7.4: Specific Speed for a Number of Hydraulic Turbines 


_ NP? _ 300 x 1000! 
S= HD/4 4005/4 = 

Comparing the values of specific speed to Table 7.2 or Figure 7.4 it is evident that a Pelton Wheel 
design is required. The exact choice of generator speed being determined by various design trade-offs 
such as between generator cost (increasing with number of poles) and windage losses (reducing with 
number of poles). 


7.4 Problems 


1. An axial flow gas turbine has three stages of compression and one stage for the turbine. The 
compressor operates with an axial inlet flow velocity of 250 m/s at 3000 rpm. The rotor 
relative exit angle is set at —34° at a mean radius of 0.75 m. Draw the velocity triangle for the 
first stage rotor exit and hence estimate the specific work output. Answer: 15.8 kJ/kg 


2. The turbine from the previous questions operates with an axial inlet velocity of 300 m/s. The 
axial velocity is kept constant through the stage. The stator blades have an absolute exit angle 
of 60° and the rotor blades have a relative exit angle of —30°. Calculate the specific work 
output and calculate the stage loading coefficient for the turbine at a mean radius of 0.75 m. 
Answers: w = 107.7 kJ/kg, Yturbine = 1.94 


3. Compare the stage loading coefficient of the turbine and the compressor in the previous two 
questions. Why is the turbine stage loading coefficient higher than that of the compressor? 
What limits compressor loading? Answer: comp = 0.28 


4. The Sellrain-Silz power station in Austria has two turbines each delivering 260 MW of power. 
The operating head is 1233 m and the rotational speed is 500 rev/min. Show that a multiple 
nozzle Pelton turbine is best suited for this application. 
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Chapter 8 


Axial Flow Machines 


8.1 Reaction for Repeating Stage 


The aim of this section is to relate reaction to blade dimensions for a common special case: a repeating 
stage with constant axial velocity. This shows how the static enthalpy drop (and therefore also the 
static pressure drop) are controlled entirely by the blade geometry. 


Consider the h-s diagram in Figure 8.1. First consider the enthalpy drop from the stage as a whole: 
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Figure 8.1: h-s diagram with ho and hore 


y2 v2 
Ahstage = hi — h3 = (roi 5 ) (ro A ) 


And this case is for a repeating stage so Vi = V3: 


Ahstage = Aho =w = Um (V26 = V30) = Um(Wa9 _ W30) 


Remember that Vg = wr + Wọ if the algebraic manipulation is not obvious. 


Now consider the rotor: Ahrotor = h2 — hg Earlier (Chapter 6) Rothalpy, J was defined: 


wW? UP? 
I = ho —UmVe = h + — - —™ 
2 2 
This was shown to be constant across a blade row, so: 
We U2 w2 
hot 2 mM i hyp 2-72 
a+ 5 5 3 + 9 Un 


Rearranging: 


1 1 1 
hə — hz = 5 (Ws w3) = 5 (Waa +V- W359 —-V;) = z Wo + Wag) (W30 — W29) 


The axial velocity V, is constant so Viz = Voz = V3x = Vz. Recall the definition of Reaction R 
and substitute for ho — hg and hı — hg to give: 


_ Ahtrotor _ haha _ (Wes + Wo9)(W39 — W20) 1 (W3 — W29) 
Alstage hy = hg Um(Wa9 = W3o) 2 Um 


R 
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We can write W39 = Vz tan 63, and Wo9g = Vz tan bə from a generic velocity triangle so: 


1 Vz 
R = “5g en + tan (33) 


Vz/Um is the flow coefficient so this expression can be simplified slightly: 


R= —5o(tan b2 + tan 63) (8.1) 


It is very useful to change Equation 8.1 into one that uses a2 instead of (2 this form of the 
equation is particularly useful as it uses the stator geometry. 


Apply the key rule in the tangential direction so: 


Wag = Vag — Um = Vz tan ag — Um 


Insert expressions for W29 and W39 into the previous equation for reaction: 


1 (W39 — W: 
R=-5 ( in 20) 
Substitute into Equation 8.1 to get: 
R= ie (tan a2 + tan 83) + ae 
which simplifies to: 
R= ; — (tan a2 + tan 33) (8.2) 


Example A turbine stage of a reaction 0.4 is to operate at a flow coefficient of 1.0, if the stator exit 
flow angle is 60°. Estimate the rotor relative inlet and exit flow angles. 


Solution This example is quite easy but drives home the point that once the non-dimensional flow 
parameters are set the geometry of the machine is fixed. 


The Stator exit flow angle is known ag along with the reaction, R and the flow coefficient ¢ 
Rearrange Equation 8.2: 


1—2R 
tan 03 = J — tan a2 
and then substitute the values we are given: 
1—2 x 0.4 


B3 = tan! ( — tan 60°) = —56.9° 


1.0 


We now examine two special cases of zero reaction (also called impulse blading) and 50% re- 
action to illustrate the fundamental trade-offs that are made when varying the level of reaction. In 
practise machines will operate with varying levels of reaction but comparing two very different cases 
gives information about the general trends. 
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8.1.1 Zero Reaction (Impulse) Stage 


If R = 0, from Equation 8.2 it is apparent that tan 32 = — tan (33 So the inlet and outlet of rotor have 
equal and opposite and opposite angles. This leads to bucket shaped blades. 


|W3| = |Wo| and 63 = -62 


A sketch of an impulse blade is shown in Figure 8.2 on the left hand side. We saw earlier that for 
a rotor row Equation 5.20 that energy conservation is equivalent to: 


U2 U2 
horelt = = = horel2 = T 


Since there are no changes in blade speed between the inlet and the outlet of the rotor then 
horel1 = Rorelg which means that for the zero reaction case: 


w2 we 
horen = hore => hi + — = h + 


7 a m a 


so there is no static enthalpy change across the rotor in an impulse blade. Since: 


oe AProtor 


Ria SS 
ADstage 


there is no pressure change across rotor and all acceleration of the flow is across the stator, which 
is why it is often called a nozzle. Note that: 


e Zero reaction has high turning in the rotor, typically 62 = 60° = — 63 giving around 120° of 
flow turning in rotor. 


e All the stagnation enthalpy drop must be across rotor, since only the rotor is moving. There is 
no work done in the stator. 


Zero reaction is often called impulse blading as there is no pressure change in the rotor so the 
mechanism of turning the rotor is down to kinetic energy changes. Where R > 0 there is a pressure 
drop in the rotor passage so the mechanism of energy transfer is more complex. The Euler equation 
however Equation 5.12 shows that the detailed mechanisms within the blade row do not determine 
the work input and output from the stage - it is the change in angular momentum through the rotor 
passage that matters. 


8.1.2 50% Reaction Stage 


It is useful to examine what happens when the reaction level is 50%. 


Equation 8.1: 
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Impulse Blading 50% Reaction Blading 


Figure 8.2: Impulse and 50% Reaction Blading 


(tan ag + tan 83) 


oe b3 = —Q2 
and since 
tan &3 = tan 3 + 2 
5s 3 V, 
and 


2 2 V 


tan @2 — tan b2 = tanag — tan b3 


=> tan. = — tanas =>> ~o =—a3 


Now a3 — q; as the stage is a repeating stage. This means that the rotor blades are mirror image of 
the stator blades. 50% reaction gives equal acceleration through rotor and stator, and equal pressure 
drop. 


The two types of blade design are sketched in Figure 8.2, this illustrates that when presented with 
a blade geometry from a simple examination of the inlet and exit angles the level of reaction chosen 
for the particular blade can be estimated. Blades with high turning are generally impulse blades. 


8.2 Loading and Efficiency Variation with Reaction 


Changing the reaction (for constant V, and a repeating stage) will change the blade geometry (Figure 
8.2). What implications does changing blade design have on loading and performance? Consider a 
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turbine stage with constant axial velocity, repeating stages and an axial exit velocity. For impulse 
blades R = 0 and we know from Equation 8.1 that the blades will be bucket shaped. 


So W2 = —W39 and from the exit velocity triangle V3g = 0 = U = —W3¢ 


So at position 2: 


V29 =U + Wo = U + U = 2U 


Euler Equation: 


w = U (Vz — V29) = U2U = 2U? 


So W = 2 for impulse blading, repeating stages, constant axial velocity and zero swirl at exit. 
For a 50% Reaction design where R = 0.5 and we have G2 = a3 = 0 and $3 = —az 


Rotor/stator velocity triangles are a mirror image so from the velocity triangle at 2: 


Və = U 
From the velocity triangle at station 3: 
V39 = 0 
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Euler Equation: 


So Y = 1 for 50% reaction blading, repeating stages, constant axial velocity and zero swirl at 
exit. So there is twice as much work for an Impulse Stage compared to a 50% Reaction Stage for the 
same design constraints. 


8.3 Stage Efficiency 


Turbine efficiency is related to the energy loss or entropy gain within turbomachines. This is an 
enormously complex topic and the best (although very advanced!) reference is Denton (1993). There 
are three basic mechanisms for entropy creation inside a blade row: 


1. Viscous friction in either boundary layers or free shear layers. 
2. Heat transfer across a finite temperature difference. 


3. Non-equilibrium processes such as occur in shock waves. 


The latter two are complex subjects and left aside in this book. The entropy rise in a boundary 
layer can be taken to be proportional to a Reynolds number based on boundary layer thickness. In 
particular for a turbulent boundary layer the dissipation coefficient (a measure of the rate of entropy 
generation in the boundary layer so smaller values are good). 


Cp = 0.0056Re, ° (8.3) 


where Re, = ere where p, u are fluid properties and V is the free-stream velocity and b is the 


boundary layer thickness. The origin of equation 8.3 is a subject for advanced fluid mechanics but 
the implications can be discussed here. For a given boundary layer thickness b a high free-stream 
velocity will result in a greater value of Re, and therefore a greater value Cp, more entropy and 
therefore lower efficiency. 


Simplifying fluid dynamics slightly we might therefore expect the stage efficiency to be somewhat 


dependent on the velocity in the machine - the higher the velocity the lower the machine efficiency. 
Consider the velocity value at stator exit V2. 


e 0% Reaction, Impulse, V2 = 4/ V2 + V? = y V2 + (2U)? 


e 50% Reaction V2 = 4/ V2 + VŽ = V V2 + U? 


So for a comparable blade the higher velocity for impulse blading gives more frictional loss and 
lowers peak efficiency. The choice of reaction therefore is fundamentally a trade-off between high 
loading and high efficiency. 
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Figure 8.3: Locations for Tip Clearance Flow 


8.4 Choice of Reaction for Turbines 


In any situation with a rotating blade and a stationary casing there has to be a small “tip clearance” 
between the end of the blade and the casing in which the machine rotates. Since the tip of the blade 
has a pressure difference over it there is leakage flow over the top of the blade, the flow here is highly 
complex. In a turbine the flow over the tip of the blades does no work so tip leakage flow is usually 
detrimental and should be minimised as far as possible. ! Tip leakage losses occur at the gap between 
the stator tip and rotating hub and at the gap between the rotor tip and casing clearance as shown in 
Figure 8.3. 


Leakage flow driven by the Ap across each row which is linked to reaction as from equation 6.3.: 


R= Ahroror _, AproTor 
AhsTAGE APsTAGE 


For 50% reaction, there is a pressure drop across both the stator and the rotor. So we get some 
leakage for each row, but the sealing is relatively straightforward as the pressure drops are quite 
modest. So we can make do with a simple “drum” construction. (Top half of Figure 8.4) 


For 0% reaction (impulse) all the pressure drop is across stator so the leakage flows would be 
much larger. A “Disc and Diaphragm” construction is often used. This enables the leakage path to 
be lowered to a small radius so although Ap is quite large the area for flow is reduced, reducing the 
leakage. However it requires a more complex construction. (Lower half of Figure 8.4) 


Gas conditions might also influence the choice of reaction. At turbine inlet, the fluid temperature 
and pressure are very high, high work per stage reduces this temperature and pressure quickly. There 
are a number of reasons this might be beneficial. High work per stage may enable the designer to 
avoid the expense of having an extra stage in the machine reducing the build cost. The materials 


'For a compressor one generally wants to minimise tip leakage but the situation is much more complex and sometimes 
a small amount of tip leakage can help the overall stability of the machine. 
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Figure 8.4: Schematics of Disc and Diaphragm Construction 


required to operate at high temperatures for both steam and gas turbines are incredibly expensive so 
reducing the machine temperature rapidly may again reduce the manufacturing costs. Finally most 
high performance gas turbines have very complex blade cooling arrangements as the temperature of 
the fluid going around the turbine blades is much higher than the melting point of the material they 
are made off - requiring only one stage of this blading has clear advantages in terms of reduced design 
effort and manufacturing cost. 


8.5 Compressor Design 


Compressors operate in an adverse pressure gradient, that is the exit pressure is much higher than the 
inlet pressure. There is a limit to how much pressure rise you can obtain in a single stage without 
boundary layer separation. Boundary layer separation means that the pressure rise through the ma- 
chine is removed and if this happens suddenly you have the temporary situation where there is a high 
pressure in the combustion chamber and no pressure rise through the compressor to prevent the flow 
of air from high pressure to low pressure. This leads occasionally to flames coming out the front of 
aircraft engines! Though modern design much reduced the incidence of these events. 


So for a turbine excursions from design values of ¢ or W are likely to reduce efficiency, for a 
compressor they are likely to compromise the operation of the whole machine. This means that 
the design imperatives for compressors are about ensuring consistent performance, i.e. avoiding 
phenomenon such as stall and surge rather than achieving maximum efficiency. For compressors the 
values of stage loading and flow coefficient are therefore more modest than for a turbine and the level 
of reaction tends to be around 50% to ensure that the adverse pressure gradient is shared between the 
rotor and the stator as far as possible. 
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8.6 Multistage Steam Turbine Example 


Here many of the techniques developed for axial flow machines are illustrated using an axial flow 
steam turbine design example, though the same technique could be employed for gas turbines or 
compressors. 


Example A steam turbine is to be supplied with superheated steam at conditions of 10 bar and 
400°C the exhaust pressure is set to 4kN/m?. The turbine isentropic efficiency is expected to 
be around 88% based on past experience with machines of this type. The machine will rotate at 
3000 rpm, a speed set by the electrical grid frequency. The mean radius of the machine will be 
0.75 m. 


Repeating stages, with constant axial velocity and axial leaving velocity are to be used. Estimate 
the number of stages and make a choice of blading with a stator exit angle of 75°. 


Solution This is a design type problem - so there is no one right or wrong answer. The first thing 
to be clear about is the difference between the changes in steam conditions over a stage and over the 
machine. The inlet conditions for the machine are fixed by the delivery pressure and temperature of 
the inlet steam. The exit pressure at the machine exit is fixed by the condenser pressure. The pressure 
drop of the machine what the blading we choose has to deliver. A good first step is to examine the 
machine performance. 


To make clearer the distinction between machine conditions and stage the entry and exit condi- 
tions for the machine are given the labels a and b respectively. The steam conditions at inlet and exit 
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are determined by reference to Steam Tables or charts (For example Rogers and Mayhew (1994). At 
inlet the steam conditions are: 


ha : 10 bar, 400°C = 3266 kJ/kg 


if we then sketch a line of constant entropy on the Mollier chart from point a to the 0.04 bar line or 
interpolate using tables we can determine hys: 


hbs = 2253 kJ/kg 


The change in enthalpy across the machine is therefore: 
Ah = ha — hos = 3266 — 2253 = 1013 kJ/kg 


Now: 
n= Ahactual 
Ahideal 


This is the enthalpy drop across the machine. The mean blade speed is given by: 


=> Ahactual = nAhideal = 0.88 x 1013 = 891 kJ/kg 


20 
Um = wrm = 3000 x 60 x 0.75 = 236 m/s 


Firstly consider a design based on impulse or zero reaction design. Earlier we saw how the stage 
loading coefficient for impulse blades with repeating stages, axial velocity and zero swirl at exit was 


p=2 


w 
UZ, 
this is the work output from each stage of the machine. The number of stages is given by the machine 
enthalpy drop divided by the stage enthalpy drop or: 


891 
— = 8.03 
111 


This would be taken as eight stages for impulse blading. 


w wU? = 2 x 2367 = 111 kJ/kg 


Nstages = 


The velocity triangles for both impulse and reaction blades are shown in Figure 8.2 on page 94. 
For a stage we have inlet designated by station 1, stator exit by station 2 and rotor exit by station 3 in 
contrast to the a, b labelling for the machine. 


From the velocity triangle for an impulse stage: 


V9 472 
Vog Um 72 m/s and V2 ano nT m/s 
Again from the velocity triangle: 
Wao V29 = Um Um = Um —1 236 Ə 
an Be Vox Vox Vox p2 = Vox = 126 


Since the blades are “bucket” shaped G3 = — 82 = —61.9°. The deflection or turning across the rotor 
is therefore: 2 x 61.9° = 124°. The key point is that once the level of reaction has been specified this 
determines the blade geometry! 


Secondly consider a design based on reaction blading where R = 0.5. The stage loading coeffi- 
cient for reaction blades with repeating stages, axial velocity and zero swirl at exit was ọ = 1 


Y 


w 
2 
Um 


w = WU? = 1 x 2367 = 55.5 kJ/kg 
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this is the work output from each stage of the machine. The number of stages is given by the machine 
enthalpy drop divided by the stage enthalpy drop or: 


891 
tages — 55.5 


This would be taken as sixteen stages for reaction blading. 


Ns = 16.05 


The determination of blade geometry for reaction blades is straightforward, a cascade view and 
velocity triangle are shown in Figure 8.2. We know that the rotor blades are the mirror image of 
the stator blades. Therefore G2 = —æ3 = 0 as there is no exit swirl. Stator and rotor blades are a 
mirror image: 33 = —@2 = —75°. The deflection across the rotor is therefore 75°. Compared to 
the impulse blading the 50% reaction blading gives lower turning and therefore lower work output, 
necessitating more stages for the same power output. The advantage of reaction blading is a slightly 
higher efficiency. 


In practise the design is not limited to a single choice of reaction, so one might end up with four 
impulse stages in a high pressure turbine and eight stages of 50% reaction in a low pressure turbine. 


Example In the previous example the turbine is designed to produce 100 MW of power. At inlet 
the steam density is around 3.26 kg/m? (this can be checked with steam tables). Estimate the local 
blade height h. 


Solution The first step is to find the mass flow. 


P 100~x 10° 
P= n=- endear ho) 
a m= w 891 x 108 oji 
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From continuity (mass conservation): 


m 
m = pV, p42, 2) SS ae 


The axial velocity V, = 126 m/s and the steam density p = 3.26 kg/m? so: 


112.2 


Ee a 


The rationale behind these examples is to show that even with a relatively small amount of theory 
one can obtained powerful insights into the reasons behind the design choices made by turbomachin- 
ery designers. 


8.7 Problems 


1. A steam power plant has a High Pressure (HP) turbine, an Intermediate Pressure (IP) turbine 
and a Low Pressure (LP) turbine. All stages of the HP, IP and LP turbines are designed for an 
axial leaving velocity and a rotational speed of 3000 RPM. The axial velocity is kept a constant 
through each turbine. 


e HP turbine: stagnation pressure drop from 15M N/m? to 3M N/m?; isentropic turbine 
(not stage) efficiency (total-to-total) of 0.88; impulse blading; stator exit flow angle of 
75°; blade mean radius of 0.51m. 


e IP turbine: stagnation pressure drop from 3M N/m? to 0.3M N/m?; turbine isentropic 
efficiency of 0.88; 50% reaction blading; stator exit flow angle of 70°; blade mean radius 
of 0.74m. 


e LP turbine: stagnation pressure drop from 0.3M N/m? to 4kN/m?; turbine isentropic 
efficiency of 0.84; 50% reaction blading; stator exit flow angle of 70°; blade mean radius 
of 1.09m. 


Superheated steam enters the HP turbine axially at a stagnation temperature of 550°C’. The 
steam is reheated to 550°C’ between the exit from the HP turbine and the inlet to the IP turbine. 


(a) Use the steam tables or chart to find relevant energy changes and estimate the number of 
stages required in the HP, IP and LP turbines (the kinetic energy leaving the LP turbine 
may be neglected). 


(b) Calculate the axial velocity for the HP, IP and LP turbines 


(c) Estimate the relative flow angles at exit from the rotor blades for each turbine 


Answers: 1) 8 stages, 11 stages, 5 stages, 2) 86 m/s, 85 m/s, 125 m/s, 3) —61.8°,—70°,—70° 


2. The repeating stage of an axial flow steam turbine is designed with the following flow angles: 
Exit flow from stator, a2 = 75°; relative flow at inlet to rotor, G2 = 45°; relative flow at outlet 
from rotor, G3 = —68°. Sketch the velocity triangles at the stator exit-rotor inlet and at the 
rotor outlet. 


3. For the machine in the previous question Show that the values of flow coefficient and stage 
loading coefficient for the stage are o = 0.366, Y = 1.272. 
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4. An intermediate pressure steam turbine is to be designed with 7 repeating stages as described 
above. The inlet conditions are: total temperature = 550°, total pressure = 4.0 MN/m?. The 
exit total pressure is 0.5 MN/m?. The rotational speed of the turbine is to be 1500 rev/min 
and the power output required is 90 MW. It is estimated that the turbine efficiency will be 
90%. Calculate the mass flow through the turbine, the axial velocity, mean blade speed and the 
mean diameter of the turbine. Answers: 162 kg/s; 91.4 m/s; 249.7 m/s; 3.18 m 
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Chapter 9 


Hydraulic Turbines 


Hydro-electric power accounts for up to 20% of the world’s electrical generation. Hydraulic turbines 
come in a variety of shapes determined by the available head and a number of sizes determined by the 
flow rate through the device. This chapter introduces no new concepts but instead allows the analysis 
techniques learnt over the preceding chapters to be applied by means of analysis examples. 


Three three basic hydraulic machines (Pelton, Francis and Kaplan) are discussed and with all 
hydraulic machines the simplification for analysis and operation is that the fluid density is constant 
but the complication is that cavitation may occur in the machine. Cavitation is effectively the boiling 
of the liquid as it is exposed to extremely low pressures in certain parts of the turbine. 
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Figure 9.1: Pelton’s Patent Application and Analysis Model 


9.1 Pelton Wheel 


This device invented around 1880 by Lester Pelton is used in high head locations with 300m to 4000m 
head, it has a peak efficiency of around 90%. 


Figure 9.1 shows the Pelton wheel on the left and the sketch on the right shows a diagram of 
operation with velocity triangles. 


The analysis of this machine is straightforward and can be found in any basic fluid mechanics 
textbook (e.g. Massey (1989)) as a simple control volume can be drawn around the rotor blade and 
the force determined from a simple analysis. Alternatively the results can be obtained from the Euler 
equation as follows. 


Recall that P = mw(Vogr2 — Vigr1) so the key parameters to work out are the inlet tangential 
velocity and the exit tangential velocity. 


At inlet the situation is straightforward Vig = V; the jet velocity, as the nozzle directs the flow 
only in the tangential direction. At exit the velocity triangle shown in Figure 4.11 must be used. 


Note that: W; = V; — U so W2 = kW, = k(V; — U) where k is an empirical coefficient for 
frictional losses in the bucket. From the velocity triangle: 


Vag = U + Wag 


and that W29 = —W2cos O, where © is the bucket angle shown in Figure 9.1. angle made by 
the bucket. The negative sign arises as our sign convention is “positive in the direction of rotation”. 
Substituting yields: 


V29 = U — W2 cos © = U — k(V; — U) cos O = —kv; cos O + U(1 + k cos ©) 
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Substitute this into the expression for power: 


P = mU (V; —U)(1+kcos@) 


In the Pelton wheel all the pressure drop occurs in the stator or nozzle so the machine can be 
classified as an impulse machine. 


9.2 Analysis Approach 


The analysis of the Pelton wheel is extremely straightforward, other hydraulic machines required a 
more nuanced approach. This is a three step process: 


1. Given the flow rate apply the continuity equation to get the radial or axial velocity. 


2. The geometry (blade angles and radii) will yield absolute and relative velocities by means of 
velocity triangles. 


3. Power is obtained from the Euler Equation. 


Empirically obtained values of efficiency (or estimated values based on prior experience) link 
the power produced and the head drop across the machine. There are two methods of doing this 
depending of what information is available. 


1. Ifthe total head is given this yields the ideal power: 
Prdeat = pQgH = thgH 


Given the ideal power, we have the actual power (from the Euler equation above) this yields 
the efficiency. 


2. If loss information is given, the actual head is obtained from: 


Pactual 
actual = Hioss 


pQg 
The ideal head is given by Hactual — Hioss so given the actual power and the losses efficiency 
can be obtained. 


This may seem slightly abstract at the moment but will become more concrete as the examples are 
explained. Firstly a Francis turbine analysis is conducted where the overall head is given and the draft 
tube and supply pipe losses must be estimated, secondly an analysis on a Kaplan turbine is presented 
where losses are calculated for each component. 


9.3 Francis Turbine 


The Francis turbine was briefly described in Chapter 4. But for this analysis of the machine each 
station through the device is considered in more detail, a diagrammatic representation of the turbine 
is shown in Figure 9.2 with each of the analysis stations labelled. The four stations are: 
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Datum height Oa 
tailrace 


Figure 9.2: Analysis of a Francis Turbine 
1. Inlet to Guide Vanes. (Station 1) There may be a small amount of swirl at inlet to the guide 
vanes because of the use of a volute to distribute the flow 


2. Exit from Guide Vanes. (Station 2) Flow is deflected by the guide vanes (stator) to give high 
swirl velocity, and thus high angular momentum. Draw the velocity triangle to get the relative 
flow at inlet to the runner blades. 


3. Exit from Runner. (Station 3) The aim here is to design at exit for zero absolute swirl as 
otherwise kinetic energy is wasted. Again we must draw the velocity triangle to get relative 
flow at exit from runner. 


4. Draft Tube. (Station 4) This diffuser allows recovery of some of the exit kinetic energy of the 
axial flow. 


Unlike the Pelton wheel all the pressure change does not occur over the stator so the Francis 
turbine can be thought of as a reaction machine. Flow from station to 1 station 2 through guide vanes 
is radially inwards, while being deflected and gaining angular momentum. 


First we apply continuity: Flow rate Q; = Qa: 


Q2 = Vər 2Tr2b2 and Qı = Vir2rrıbı 


If bı = be, (constant width vanes), then Vərr2 = Virrı. Now the guide vanes have a non-zero 
thickness, and this introduces a blockage t at station 2 so to allow for this: 


Q2 = 2mrebe(1 — t) Vər 


Typical t = 0.08 or about 8% of the area is blocked by the inlet guide vanes. The detailed analysis 
process is as follows: 
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e From the value of Vər, the flow angle at exit from the vanes, œz and the blade speed, wRa, 
construct the velocity triangle at vane exit, 2 


e Usually the design is for no exit swirl to reduce the exit kinetic energy that is lost, i.e. V39 = 0. 
Thus from the blade speed wRg, the velocity triangle at runner exit, 3 can be constructed. 


e The runner inlet and exit relative flow angles (to give the geometry of the runner blades) are 
obtained and also the power from the change in angular momentum across the runner (Euler 


equation). 


e Finally to get the efficiency, the height of the reservoir and information on losses is required. 


Example Consider a machine with the following specification: 


e Outer diameter of runner, 2r2 = 2.0 m 

e Rotational speed, N = 200 rev/min 

e Guide vane height, b2 = 0.3 m 

e Vane blockage, t = 0.08 

e Vane exit angle, œ = 75° from the radial direction 

e Impeller designed for axial flow at exit, ag = 0, and r3m = 0.5 m; b3 = 0.4m 

e Supply head, Ho = 63 m; Flow rate Q = 12 m?/s 

e Flow losses: 2 m of heat loss in the supply pipe, 0.5 m of head loss in the draft tube 
e Draft tube velocity of 4 m/s 


Determine the relative angles at inlet and exit of runner to give a preliminary design for the runner 
geometry. Also find the power output. 


Solution The strategy is to apply the analysis steps outlined earlier. The most difficult part is 
drawing the velocity triangles at station 2 and station 3. 


Continuity at 2 At exit from guide vanes: 


Q = 2nrgbo(1 — t)Vo, => 12 = 2r x 1.0 x 0.3(1 — 0.08) V;2 


Working the numbers though yields: Vər = 6.92 m/s 
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VAA N 


Figure 9.3: Velocity Triangle for Francis Turbine Guide Vane Exit 


Velocity Triangle at 2: Figure 9.3 shows the velocity triangle at 2 and the corresponding cascade 
view. Earlier velocity triangles were drawn for pumps in a very similar manner. Since the geometry 
of the blades is very complex only the first part of the blade is drawn. To determine power the Euler 
equation only requires values at inlet and exit of the turbine so a sketch at inlet and exit will suffice. In 
reality the turbine blade curves around into the axial direction which can be represented using CAD 
for example but is not easy in a sketch. The velocity triangle is in the tangential radial plane. 
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Vo9 = Vər tan az = 6.92tan 75° = > Vog = 25.83 m/s 


Vo = 1/ V2 + V2, = v 25.83? + 6.922 = 26.74 m/s 
2r 2 


From the velocity triangle: 
V29 = Wag + wre 


200 
wre = 2T X a- 1.0 = 20.94 m/s 


=> Wag = Vag — wrz = 25.83 — 20.94 = 4.89 m/s 


So we can now work out the relative flow angle at 2: 


W. 4.89 
tan By = —* Bo = tan`! (= =) = 35.2° 


In Chapter 4 the velocity of a jet expanding to atmospheric pressure was shown to be \/2gH for 
this turbine the ratio, V2/\/2gH = 0.76 so not all the acceleration of the fluid occurs over the stator 
so the reaction R > 0 and the machine is not an impulse turbine. 


Continuity at 3: The flow area is given by 277T3,,,b3 so: 
Q = 277 3mb3V3x 
Rearrange and put the numbers in to give V3,: 


V3q = 12/ (20 x 0.5 x 0.4) = 9.55 m/s 


Velocity Triangle at 3: Figure 9.4 shows the velocity triangle at 3 and the corresponding cascade 
view, again the complex geometry in the radial direction is not drawn instead only the exit of the 
runner is considered. 


No exit swirl, V39 = 0 


200 
<. W39 = —wr3m = —27 X a x 0.5 = —10.47 m/s 


So we can now work out the relative flow angle: 


Power output: Recall that: 
P = thw (r2 V29 — T3mVo3) 


but V39 = 0 so: 


200 
P = mwr2V2 = 1000 x 12 x 2r x G0 x 25.82 = 6.49 MW 


114 


Download free eBooks at bookboon.com 


Basic Concepts in Turbomachinery Hydraulic Turbines 


X 


VAAN 


3 
© PIN Ww, 


runner exit 


Figure 9.4: Velocity Triangle for Francis Runner Exit 


Draft Tube Analysis Although external to the rotating parts of the turbine the draft tube is an 
important part of the hydraulic machine. The draft tube is a conical diffuser with around 7° divergence 
which reduces the exit kinetic energy in the departing fluid and therefore increases the efficiency of 
the machine as a whole. 


In this example the aim is to reduce exit kinetic energy at station 4, note that h4 is normally below 
the height of the river that the flow exits into. 


As seen in Chapter 4 the total head at station 0 (the reservoir) is: 


Po , Vo 


+z 
pg 2g °° 


But po = 0 as in hydraulics we always use gauge pressure where zero is the ambient atmospheric 
pressure, in the reservoir Vo = 0 so the total head at station zero is simply given by Ho = zo 


We need to find the head across turbine: 


AH = H; — H; 


Now Hı = Ho — h fp, where h pp is the pipe head loss due to friction. Now examine the total head 
from the other end of the machine starting at station 4. If loss in draft tube is given by h fpr then: 


H, = H3 — hfjpr => H3 = H4 +hfpr 


The total head at station 4 is given by: 


V2 
mang ig 
pg 2g 
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The pressure at station 4 must be equal to the pressure of the fluid in the river as a whole so is 
given by the pressure due to a height of fluid: p4 = pgh4. Also z4 = —h4 substituting gives: 


h VŽ VŽ 
— PI = 4: 4 h4 = 4 


H4 
pg 2g 2g 


Thus 
2 


Beime 
3 fDT 2g 


So head across turbine: 


V2 
Hı — H; = Ho — h fp — (ror + x) 


simplifying slightly: 
2 


V, 
Hı — Ha = Ho- 35 — hyp — hype (9.1) 


So for maximum head change across turbine and therefore the maximum power the design aim is 
to keep V4, h fp and h fpr as low as possible. 


To apply this information on pipe flow losses and draft tube losses is required and the velocity at 
station 4. In this example the information is provided: hs, = 2 m, hepr = 0.5 m, Vy = 4m /s. 


h fp can be calculated from a simple knowledge of pipe flow friction (Massey (1989)). Draft tube 
losses require some empirical knowledge and V4 is set by conservation of mass. In this example 
V3 = V3z = 9.55 m/s, so the exit velocity is considerably reduced by the draft tube. 


Apply Equation 9.1 with the above numbers: 


2 


4 
H= = Od — 5 Am = Gam 
g 


The turbine efficiency may then be worked out: 


Pactual 6.49 x 10° 
Pyaea) 1000 x 12 x 59.68 x 9.81 


9.4 Kaplan Turbine 


The Kaplan turbine is shown in diagrammatic form in Figure 9.5 which also shows the analysis 
stations. In comparison to the analysis of the Francis turbine an extra analysis station is needed 
(station 3) between the guide vane exit and the runner inlet as the flow is turned through 90° before 
entering the turbine blades. The analysis approach is as follows: 


1. From station 2, the guide vane exit, to station 3, the runner inlet one can either use conservation 
of angular momentum if guide vanes are radial. In some Kaplan turbines the guide vanes are 
orientated in the axial direction in which case station 2 = station 3. The velocity triangles are 
then drawn to get relative velocities at rotor inlet. 
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Figure 9.5: Analysis of a Kaplan Turbine 


2. At 4, runner exit, draw the velocity triangle to obtain relative and absolute exit velocities. 
3. The Euler equation 5.12 gives the power P = mw (r3 V303 — r4 Vag) 


4. Calculate losses. In this example component losses will be considered. Component head losses 
will be obtained by working through the stations from 0 to 4. This gives the turbine efficiency. 


5. Finally the draft tube, station 4 to station 5 is analysed, this calculation will give us the overall 
efficiency of the site. 


Example Consider a Kaplan Turbine with the following characteristics: 
e Guide Vane Exit: Radial flow, r2 = 2.0 m; blade height, bo = 1.0 m, blockage, tə = 0.08; 
absolute flow angle to radial, a2 = 50° 


e Runner: Mean radii at station 3 and 4 are the same: r3m = Tam = 0.85 m, blade height, b3 = 
0.7 m; designed for zero exit absolute swirl, V4g = 0; rotational speed, N = 300 rev/min. 


e Draft Tube: Reduces axial velocity to V; = 5 m/s 


e Flow rate: Q = 36 m3/s 


Calculate the runner relative and absolute flow angles and velocities at inlet and exit. Also calcu- 
late the power output. 


Solution The strategy is to apply the analysis steps outlined earlier. In a Kaplan turbine the flow 
through the guide vanes may be be radial or axial so the geometry must be carefully inspected. 
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Figure 9.6: Velocity Triangle for a Kaplan Turbine at Guide Vane Exit 


Guide Vane Exit to Runner Inlet At guide vane exit (station 2), the flow is radial. Apply 
continuity: Q = 2mrgb2(1 — t2)Vər, where Vər is the radial velocity at 2 so: 


36 


Vor = 
2r — Om x 2.0 x 1.0 x (1 — 0.08) 


= 3.11 m/s 


From the velocity triangle at 2 (Figure 9.6): 


V29 = Vər tan ag = 3.11 x tan 50° = 3.71 m/s 
V2 = (Vor)? + (Vag)? = y 3.112 + 3.712 = 4.84 m/s 


Since there are no blades between station 2 and station 3, angular momentum is conserved, so: 


r2 V: 2.0 x 3.71 
T3m V39 = r2 V29 V30 - —— <-> 8.73 m/s 
i 


At station 3 the flow is axial. The continuity equation is: Q = 277r3mb3V3z, where V3z is the 


axial velocity at 3. 
36 


Ve, = 
3e On x 0.85 x 0.7 


= 9.63 m/s 


From the velocity triangle at 3 (Figure 9.7): 
= y (Vaz)? + (V30)? = V 9.632 + 8.732 = 13.00m/s 
V: 8.73 
a3 = tan! (2) = tan! Ga = 42.2° 


The blade blade speed U3 is given by: 


27 N 27 x 300 
m — 
60 3” 60 


WT3m = 0.85 = 26.70 m/s 
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Figure 9.7: Velocity Triangle for a Kaplan Runner 
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From the velocity triangles (or from the key rule (Equation 2.1) in the tangential direction): 
V39 =U + W39 => Wo = V39 — wr3m = 8.73 — 26.70 = —17.97 m/s 


we now have all the ingredients to work out the relative flow angle at inlet to the runner: 
_, { W30 ai /—17.97 
= tan! =t = —61.8° 
pa van ( Vas ) oe 9.63 


Across Runner At runner exit, Vig = 0. Vag = V3xz aS T4m = 13m, hence from the velocity 
triangle at 4: 


Wago = -—WT4m = —WT3m = —27.60 m/s 


1 (Wao 4 (—27.60 
= tan! =t = —70.2° 
ba an ( Vis ) an ( 9.63 ) 70 


Wa = \/ W2, + W2, = V17.972 + 27.602 = 28.38m/s 


Euler Equation 
P = thw (13V363 — raV40) 


But Vig = 0 and the fluid is water with density, p = 1000kg/m? 


P = 36 x 1000 x 26.70 x 8.73 = 8391kW 


This is the power produced by the turbine. 


9.4.1 Loss Estimation 


One way losses may be calculated is to assume that across any component the total head loss is 
proportional to the velocity head at exit: 
V2 
an=r( $) 
2g exit 


This is a concept borrowed from minor losses in pipe flow calculations, this is still an empirical 
method as one has to specify an appropriate value of k for each component. For runners the relative 
dynamic head W?/2g is used as the and the cause of most fluid friction is the the relative motion 
between the runner surfaces and the mainstream flow. 


Example With the Kaplan turbine from the previous example the loss coefficients may be taken as 
k = 0.05 for the guide vanes and bend and k = 0.06 for the runner. Ignoring pipe friction estimate 
the turbine efficiency. 
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Solution Each station is considered in turn and an appropriate value of k specified. 


e ( = 1, for this example ignore the pipe loss. Though it is relatively easy to calculate. 
e 1 = 2, for the guide vanes, k = 0.05: AH = 0.05(V7) /2g = 0.05(4.847)/2g = 0.06 m 


e 2 => 3, k = 0.05 for the bend from radial to axial flow: AH = 0.05(V?)/2g = 0.05(13.07)/2g = 
0.43 m 


e 3 = 4, for the runner use the relative exit head. With k = 0.06: AH = 0.06(W?)/2g = 
0.06(28.387)/2g = 2.46 m 


From station | to station 4, across the whole turbine, the ideal total head drop for the actual power 
is given by: 


j P 
mmgAHideal =P ie) = AHi = actual =. 
mgA Fideal 
this means that: 
8391 x 10 
AHideal = = 23.76 m 


36 x 1000 x 9.81 


which is the total head required if there were no losses in the turbine. The actual total head drop 
across the turbine, 


Hı = H4 = AHideal T 5 AHiosses 


Each of the losses have already been calculated so substituting: 


H1 — H4 = 23.76 + (0.06 + 0.43 + 2.46) = 26.71 m 


This is the required total head drop across the turbine alone to produce the desired flow rate. The 
turbine efficiency is: 


Actual Power 8391 x 108 
nT = — = = 0.89 
mg(H,—H4) 36 x 1000 x 9.81 x 26.71 


This is the efficiency for the turbine alone. 


9.4.2 Draft Tube Analysis 


Normally in hydraulic machines two efficiencies may be quoted. The first known as “turbine effi- 
ciency” excludes draft tube and inlet pipe losses from the definition of available head and gives an 
indication as to the “goodness” of the turbine runner. The second known as “system efficiency” in- 
cludes the losses in inlet pipes and draft tube and is of crucial interest to those designing or operating 
particular hydro-electric power plants. 


Example Using the Kaplan turbine in the previous example. Estimate the system efficiency if the 
draft tube loss coefficient is given by k = 0.2. 
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Solution From station 4 to 5 the draft tube reduces the axial velocity and hence the wasted exit 
kinetic energy. The total head at 5 is given by: 


V2 

Hs = 4 2425 
pg 2g 

But z5 = —hs and so ps = pghs as the pressure at station 5 must be equal to the pressure in the 
river so: 
V2 Ve 52 
Hc ha 3 + 25 = 5 = = 1.27m 
pg 2g 2g 2x 9.81 


Now AHy_s = H4 — Hs,so H; = Hs + AHy_s and k = 0.2 so the draft tube loss becomes: 


VŽ — V2 9.63? — 5? 
AH4-5 = k| =— 5 | =0.2 TA E 
4—5 ( 2g ) 0 “(S>) 0.69 m 


This gives H4 = 1.27 + 0.69 = 1.96 m so the inlet head to the turbine Hı = (Hı — H4) + H4 = 
26.71 + 1.96 = 28.67 m 


This is the inlet total head required for the turbine to deliver the actual power. So the overall 
overall efficiency for the system is given by: 


_ Actual Power _ 8391 x 10° =s 
oS mgH 36x 1000x 9.81 x 28.67 ` 
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9.4.3 Effect of Draft Tube 


To make clearer the influence of the draft tube, a calculation is carried out to illustrate what would 
happen if it was replaced with a parallel pipe, Vs = V4 = 9.63 m/s and assuming that there is no 
head loss associated with the parallel pipe. Consider what happens to the head and efficiency when 
there is a parallel pipe. 


_ V2 963 


Hy =H, = = = ——__ 
4 5 2g 2x981 


= 4.73m 


This makes the inlet head to the turbine Hı = (Hı — H4) + H4 = 26.71 + 4.73 = 31.44 m. For 
the same power output the inlet head would need to be raised from from 28.67 m to 31.43 m. The 
changed overall efficiency for the turbine is given by: 


Actual Power 8391 x 108 
ns = : = = 0.76 
mgH 36 x 1000 x 9.81 x 31.44 


So there is a substantial reduction in efficiency, the difference in cost between a parallel pipe and 
a diffuser is negligible so a draft tube is well worth having! 


There is a drawback to having a draft tube, consider the turbine exit pressure for both the parallel 
pipe and the draft tube. The total head at station 4 is given by: 


V2 V2 
Hy= y 4 424 = hat tA 
pg 2g 2g 


Where h4 is the static head at station 4: 


V2 
hg = H- + - x 
2g 


For the sake of example say that the turbine exit is 1 m above the tail race, that is z4 = 1.0m 
With the draft tube installed. H4 = 1.96 m and V4 = 9.63 m/s so the static pressure head at station 
4 is: 


2 
h4 = 1.96 2. 1.0 = —3.77 m 
2 x 9.81 


The negative sign indicates that the pressure is below atmospheric pressure. Expressing this 
pressure in Pa rather than m of fluid yields: 


pa = hapg = —3.77 x 1000 x 9.81 = —36.9k Pa ~ —0.4 bar 


Without the draft tube so with a cylindrical pipe at exit. H4 = 4.72 m and V4 = 9.63 m/s so the 
static pressure head at station 4 is: 


9.63? 
ha = 4.72 1.0 = -—1. 
Ree Surg gi oe 
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Again this is below atmospheric pressure, converting to pressure rather than head gives: 


pa = hapg = —1.00 x 1000 x 9.81 = —9.8kPa ~ —0.1 bar 


which is considerably greater. 


The problem with excessively low pressures is that they can lead to a phenomenon known as 
cavitation. Cavitation occurs if the absolute pressure falls below the saturated vapour pressure, py, 
for the water temperature. Or in other words when the pressure drops too low the water starts to boil 
even at the ambient temperature encountered by hydraulic machines. Essentially bubbles of vapour 
will be formed. As the pressure rises again the bubbles collapse suddenly and very high instantaneous 
pressures can be created (p > 500bar) which can cause significant erosion of blade surfaces and lower 
machine performance. 


In practise dissolved air starts to come out of solution of the fluid first and so some margin between 
the saturated vapour pressure and the pressure the turbine operates at is required. This is normally 
done by using an empirical parameter known as Thoma’s parameter: 


(9.2) 


where p is the pressure of the fluid, p, is the saturated vapour pressure and AH is the head drop 
across the machine. Charts of Thoma’s parameter derived from practical experience are available. 
The problem of cavitation can be alleviated somewhat by reducing z4, that is lowering the turbine 
closer to the river level or even below it. 


9.5 Problems 


1. A Francis turbine operates with the following data: Head: 100m; Flow rate 36 m°/s; Speed: 
410 rev/min; Outer radius of runner 0.6 m; Absolute flow angle: 72°; Velocity: 36 m/s; Exit 
swirl angle: 0° 


Calculate the effective runner width (i.e. ignoring blockage) and relative flow angle at inlet, the 
power output and efficiency. Answers: 0.859 m, 37.3°, 31.8 MW,90.0% 


2. A water turbine develops 480 kW with a head of 12.2 m with an overall efficiency of 86% 
when fitted with a cylindrical draft tube of 1.5m diameter. If a conical draft tube of 2.25 m exit 
diameter is substituted, calculate the increase in efficiency and power. Assume the flow rate, 
head, speed and losses in the machine are unaltered, and that there is no exit swirl. Answers: 
88.3%, 493 kW 


3. A vertical axis Kaplan Turbine is to be designed for an overall head of 9 m, delivering 800 kW 
at a rotational speed of 150 rev/min. The vanes and runner are to be mounted in an annular 
duct of outer diameter 2 m and inner diameter 0.8 m. There is to be no exit swirl, and a draft 
tube reduces the axial velocity by a factor of two. Making a reasonable assumption for the 
overall efficiency, estimate the required flow rate and hence the flow angle at exit from the 
vanes and the relative flow angles at inlet and exit for the runner. The losses of total head are 
estimated as: 


(a) across the guide vanes, 6% of exit velocity head 
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(b) across the runner, 8% of the exit relative velocity head 
(c) in the draft tube, 20% of the change in the velocity head. 


Calculate the overall efficiency (for comparison with your estimate). Answers: 10.3 m3/s, 
60.5°, —46.5°, —70.5°, 88.3% 
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Chapter 10 


Analysis of Pumps 


This bulk of this chapter builds on the introduction to pumps in Chapter 4 by providing an example of 
pump analysis. There are only a few new concepts in this chapter instead the techniques of drawing 
velocity triangles, forming the cascade view and applying the basic equations are applied to pumps. 
Before the main example the Euler equation is applied to predict pump performance, pump diffuser 
performance discussed and pump loss coefficients are detailed. 
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10.0.1 Pump Geometry and Performance 


A general arrangement diagram of a pump is in Figure 4.1 on page 44. If station 1 is inlet to rotor and 
station 2 the outlet of the rotor, then the Euler equation gives: 


P = thw(r2Vo9 — rı Vig) 


Consider the efficiency of a pump: 


Ideal Power mgH 


e = ActualPower — mw (r2Vo9 — rı Vio) 


This means that the Head H developed over the turbine is: 
H= B io(raVa =r Vio) (10.1) 


So for a large head rise one can make r2 > rı as well as relying on the change in tangential velocity 
Vo. Recall that the key rule (Equation 2.1) can also be applied in the tangential direction only: 


V29 = Wag + wre = Vor tan b2 + wre 


The relationship W29 = Vər tan G2 can be seen in Figure 4.5 or worked out from a new velocity 
triangle. The velocity V2, is obtained by applying the continuity equation at station 2: 


— l _ 
2Tr2bə(1 = t2) 


Substitute this into the expression for V29 obtained earlier to give: 


_ Q tan b2 
zZ 2Tr2b2(1 = t2) 


Q = Qrrgbo(1 = ta) Vor => Vor = 


V29 F WT2 


If there is no swirl at entry then Vig = 0. (A reasonable assumption as the inlet is usually the flow 
from a pipe and pipe flow generally has no rotational motion). Substitute expressions for Vig and Vog 
into Equation 10.1 to give: 


Np w tan b2 2.2 
H= 10.2 
g (ona go) ma 
The pump power is obtained after a short manipulation: 
P = p | —— 10.3 
P ee +w r3Q ( ) 


Equations 10.2 and 10.3 provide powerful insights into the choices designers have in determining 
pump geometry and operation. Consider Equation 10.2. The pump head (or pressure) rise produced 
increases with increasing efficiency. At zero flow rate the head produced by the pump is determined 
by the rotational speed, radius and efficiency only, because Equation 10.2 with Q = 0 becomes: 


H = z (w?r3) 
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B, +ve forward leaning 


b, = 0 Radial Exit 


b, -ve backward leaning 


Figure 10.1: Three Blade Angles at Impeller Exit 


— Constant efficiency 
H Fag Variable efficiency 


B, +ve forward leaning 


B, = 0 Radial Exit 


B, -ve backward leaning 


Q design Q 


Figure 10.2: H vs Q for Three Blade Angles 


The tan ə term is also important. Consider three cases of radial pump design: firstly where (2 is 
positive, this is often described as forward leaning, secondly where (32 is zero called radial exit and 
finally Go is negative known as backward leaning. These three geometries are shown in Figure 10.1. 


Equation 10.2 can be used to predict the performance of the pump as the flow rate varies at 
constant rotational speed. (This is generally achieved by adjusting valves at inlet and exit to the 


pump, though exactly how this is achieved is not important for this discussion). 


The three cases are sketched on Figure 10.2 which shows two lines for each of the three pump 
geometries. The first of which is with constant efficiency (blue lines) so are straight lines on the graph 
as from Equation 10.2 H is linear with Q when everything else is constant. The three cases are: 
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P 


B, +ve forward leaning 


b, = 0 Radial Exit 


b, -ve backward leaning 


Figure 10.3: P vs Q for Three Blade Angles 


e For 32 = 0 (radial exit) then the head produced does not vary with flow rate, from Equation 
10.2: tan 62 = tan0 = 0. 


e For 62 > 0 (forward leaning) the head produced increases with increasing flow rate, as tan G2 > 
0. 


e For 62 < 0 (backward leaning) the head produced decreases with increasing flow rate. 


The plots are then repeated with an efficiency that varies with the maximum value being at a 
design flow rate Qaesign and reducing efficiency either side of the design value. These are the red 
lines shown in Figure 10.2 and show a more realistic characteristic curve for the pump. The key point 
is that with a simple equation the performance of pump as a whole has been characterised. 


The power variation with flow rate can also be sketched using Equation 10.3 which is shown 
in Figure 10.3. The Q? term in equation 10.3 is significant. When (3) > 0 this term results in the 
gradient of the P vs Q line increasing with increasing flow rate. For radial exit tan 32 = 0 so the 
gradient is constant and for backwards leaning tan 32 < 0 so the gradient reduces with increasing 
flow rate. 


As a result of this backward leaning is generally preferred to forward leaning, despite the lower 
head produced. Backward leaning pumps have a higher efficiency due to lower exit swirl and forward 
leaning blades have an increasingly steep power curve. If Q is underestimated, the power requirement 
will be greatly underestimated so backward leaning blades are more tolerant of flow rate errors. 


10.1 Pump Diffuser Analysis 


For a machine with radial exit flow from the rotor, the kinetic energy may be reduced either by means 
of a vaned diffuser, or a vaneless diffuser. For a vaneless diffuser, as the flow goes out radially by 
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continuity: 
V, X r= constant = Ky or V, = Kı/r 
By conservation of angular momentum: 


Ver = constant = Kə or Vg = Kə/r 


The exit kinetic energy is therefore given by: 


V? V+V? 1 
20 2 Or 


(Ki + K3) 


So the exit kinetic energy reduces with the inverse square of the radius: 


In a flow where the total pressure is conserved there will be a a corresponding increase in static 
pressure. If this pressure rise is large enough to cause separation then vanes may be required to control 
the flow. 


10.2 Pump Losses 


Losses are often expressed in terms of fractions of a dynamic head, in a similar manner to the tech- 
nique used in the Kaplan turbine in Chapter 9. For an impeller the relative dynamic head at inlet 
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is used. This is because fluid friction causes most of the loss and the relative motion between the 
impeller surface and the mainstream flow is the cause of much of that fluid friction. The head loss for 
an impeller is given as: 


w2 
AH Imp = Kimp zy 


K mp is the impeller loss coefficient. For the diffuser the loss coefficient is based on the absolute 
exit velocity: 


vV? 
Apis = Kbit ag 


In actual fact the diffuser is the major source of loss in the pump. The loss in a diffuser may also 
be expressed in terms of a diffuser efficiency: 


A Pactual 


10.4 
APideal ( ) 


Diff = 
where Apjdeqi is the ideal pressure rise with diffusion to zero velocity, so Apjdeqi is effectively the 
dynamic pressure at impeller exit or: 


1 2 APactual P3 — p2 
APideal = V2 and npiff = I pV2 = l pV2 


where station 3 is the diffuser exit. The diffuser head loss is given by: 


AF pit ¢ = Hə — Az 


Total head is made up of pressure, velocity and height terms, assuming that the height terms are 
negligible the diffuser head loss can be rewritten: 


V; V: = V2 — V2 
AH pits = Ho — Hg = @ ie P3 3 _ P27 P3 , 2 3 
pg 2g pg 2g pg 2g 


In the diffuser V3 will be very small so assuming V3 ~ 0: 


p-p V? 
AH pits = A a 


v2 
=> p3 — p2 = pg (Z - Att ) 
pg 2g 2g cA 


P3 -P2 — P9 (2 AĦniss) _,_ AHDijf 
= ; = 
30V? oV? \ 2g VŽ /2g 
which can be written much more simply as: 
npişs = 1— Kpiff (10.5) 


which is useful if diffuser performance is expressed as an efficiency and a coefficient is needed or the 
other way around. 
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rotation 


U=or,, 


Figure 10.4: Inlet to Pump Impeller 
10.3 Centrifugal Pump Example 
Example Consider a centrifugal pump with data: 


e At impeller inlet, mean radius, rım = 0.1 m, blade height, bı = 0.1 m 


e At impeller exit, radius, ro = 0.2 m, blade height, b2 = 0.03 m, blade blockage, t = 0.08 


Backward leaning blades with G2 = —20° 


Flow rate of water, Q = 0.2 m?/s, Rotational speed, N = 1450 rev/min 


e Assume for losses: Impeller, Kimp = 10%; Diffuser efficiency = 60% 


Calculate: The relative flow at impeller inlet, absolute flow at impeller exit, the power required, 
head produced and the pump efficiency. 


Solution The basic strategy is to work form pump inlet to exit drawing the cascade view and 
velocity triangles through the machine. 


Impeller Inlet First sketch pump and velocity triangle at inlet (See Figure 10.4). 
Continuity: 


Q 0.2 
Qnrimbi 2r x 0.1 x 0.1 


Q = 2n7T1mbiVe1 => Voi = = 3.183 m/s 


Angular velocity: 
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a 


rotation 


Figure 10.5: Exit from Pump Impeller 


_ 27 N o 271450 
ain 60 


= 151.8rad/s 


Blade Speed: 


Ui = Wrim = 15.18 m/s 


From velocity triangle: 


W? = VŽ, + (Wrim)? > Wi = V3.1832 + 15.182 = 15.51 m/s 


=i 18 —1 lm —1 o 
= — Z= = T ENE = —78&.1 
By tan ( A ) tan ( A ) tan ( 1 5 8 


This defines the flow at inlet completely. 


Impeller Exit Sketch pump configuration and velocity triangles at exit, these are shown in Figure 
10.5. Continuity is applied once again: 


0.2 
= Onrebs(l — h) Vie S Va = 2597 
Q = 2nr2bo( 2) Vr2 2 Or x 0.2 x 0.03 x (1 — 0.08) mle 


Blade Speed: 


U> = wre = 151.8 x 0.2 = 30.36m/s 
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Backward facing blades with 32 = —20° and from the velocity triangle Uz + Wo9 = V2ọ but 
Wag = Wor tan Bo SO: 


V29 = Uz + Wo, tan Bo 


Now since there is no relative motion in the radial direction in the cascade view W2, = V2, and 
so: 


V29 = 30.36 + 5.77 tan —20° = 28.26 m/s 


V2 + V2 = vV 28.262 + 5.772 = 28.84 m/s 


a 20 =] 28.26 
= tan! ( 2 ) =t 23:29 \ L 78.4? 
Q2 an ( 2) an ( 5.77 ) 


which defines completely the exit flow from the impeller. 


Specific Work Input 
w = w(r2 V29 — rım V10) 


But Vg; = 0 so: 


w = 151.8 x (0.20 x 28.26) = 857.97 J/kg 


Power input is obtain directly from the specific work input: 
P = pQw = 1000 x 0.2 x 857.97 = 171.59 kW 


Head and Pump Efficiency Idea head rise is obtained from w = gAH;, so: 


Re w 857.97 


= 87.46 
g 9.81 u 


Actual Head Rise is given by AH, = AH; — losses where 


losses = AHimp + AH giz ¢ 


w? 15.51? 
Kap = 1 — nate = 1 — 0.6 = 0.4 
vV? 28.84? 


Therefore AHq;¢ = 16.96 m - note the large diffuser loss in comparison with the impeller loss. 
The total head rise across the pump is therefore: 
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rotation 


Pump inlet: 0 at atmospheric pressure 
Impeller inlet: 1 


Figure 10.6: Pump Inlet 


Ad, = 87.46 — 1.23 — 16.96 = 69.27 m 


Pump efficiency can also be calculated easily:- 


Ideal Work w; for actual AH, 
Actual Work w 


Np = 


Now w; = gA Ha minimum work required and w = gAH;. This can be confusing but recall that 
the actual work has to supply the actual head and the losses, the ideal head is therefore always higher 
than the actual head produced. The efficiency is therefore given by: 


_ gAH, _ 69.27 


m = = 0769 
e= OAH, 8746 


10.4 Net Positive Suction Head (NPSH) 


Consider the pump shown in Figure 10.6, the pump inlet station 0 is at atmospheric pressure. The 
impeller inlet station 1 must be at a lower pressure pı to ensure flow into the pump. 


For cavitation to be avoided (See Chapter 9 for details about the phenomenon of cavitation) then 
pı > Py where pẹ, is the saturated vapour pressure of the liquid in question. This is a property of the 
fluid and for water can be found in steam tables such as Rogers and Mayhew (1994). Normally this 
pressure is converted to a head h = - as pressure changes can be directly related to the physical 
location of the machine above or below the datum. 


For satisfactory operation some margin is needed between the saturated pressure vapour and the 
pressure at inlet to the impeller this is called the Net Positive Suction Head or NPSH for short: 


PTP FNP 


pg 
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Consider the total head change from pump inlet 0 at pressure pa and height z = zo = 0 to impeller 
inlet 1 where z = z1: 


Ho = Hı +hy 


where hẹ is the frictional head loss between station 0 and station 1. At station 0 the velocity 
can be regarded as very small as the reservoir that the pump is taking fluid from is generally large 
compared to the pipe inlet diameter and since zọ = 0 m the total head is simply equal to the pressure 
head which is given by atmospheric pressure so: 
2 2 
Pa pes a 
pg pg 2g 2g 


y2 
hy = ha = 2, = hş 
but NPSH = hı — hy so: 
V2 
NPSH = ha — hy F z—hy (10.6) 


This is a very important parameter in pump selection and analysis. Note how the NPSH is reduced by 
high inlet velocity, high height above the water surface and large friction losses in the inlet piping. In 
fact most pumps do not operate in the manner shown in Figure 10.6 and are usually operated below 
the water level, this means that z, would be negative increasing the NPSH. 
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In practise N PSH values are obtained using empirical criteria. One simple treatment is to use 


Thoma’s parameter given by: 
NPSH 
g = 


AH 


(10.7) 


where AH is the total head rise across the pump, cavitation is deemed to occur when o = Cerit 
which has to obtained empirically. 


10.4.1 Cavitation Example 


Example Consider the previous centrifugal pump example. Calculate the impeller inlet height 
above the free surface water level for no cavitation. If the critical value of Thoma’s parameter is 
Ocrit = 0.1. The water temperature is 25° and atmospheric pressure = 1.0 bar. The inlet duct loss is 
25% of the inlet dynamic head to the pump. 


Solution To determine the maximum applicable height the inlet loss is examined: 


2 


v? vi 
AH; = ho- | h + 5> +a) => za = ho- h -> — AH; 
2g 2g 


The problem then reduces to finding values for ho, h1, Vi and Ah;. The NPSH can be calculated as 
follows: 


NPSH = hy — hy = Cerit x AH = 0.1 x 69.27 = 6.93 m 


Now py = 3.17k Pa at 25° for water, this is obtained from standard steam tables (Rogers and Mayhew 
(1994)). This pressure is converted to a head as follows: 
3.17 x 108 


Pv 
A = 0.323 m(ab 
” = og 1000 x 9.81 riLaba) 


Absolute pressures and not gauge (relative pressures) must be used. From the definition of NPSH: 


hı = NPSH + h, = 6.93 + 0.323 = 7.25 m(abs) 


The inlet loss is 25% of VÊ /2g so: 


VA 3.183? 
1 = 0,25 


AH; = 0.25 — 
i 2g 2 x 9.81 


= 0.129 m 


At the inlet to the pump (station 0) the total head is given by the inlet pressure head as V3, zo = 0 
SO: 
Da 1.0 x 10° 


h — — 
° pg 1000 x 9.81 


= 10.194 m(abs) 


So the maximum height above the water at which the pump can be placed is given by: 


2 


V 
zı = ho — hı — — AH; = 10.194 — 7.25 — 0.516 — 0.129 = 2.299m 
g 


So the maximum height is around 2m above the water. 
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10.5 Application to Real Pumps 


This book is about the basic concepts in turbomachinery and if you were to design a real pump you 
would need to take account of more advanced concepts such as slip - where the exit from the pump 
blades does not exactly follow the angle set by the blades. Although there is not space for it in this 
book for further details see Dixon (2005). 


10.6 Problems 


1. A radial flow pump delivers 0.015 m/s when running at a speed of 1470 rev/min. The 
impeller has a uniform width of 15 mm, with an entry diameter of 70 mm and an exit diameter 
of 120mm. It has backswept blades with a relative exit flow angle of —50°. The blade blockage 
at exit is 8%. Calculate the power required to drive the pump and the relative flow angle at inlet 
to the impeller. Answers: 803 W;—49.8° 


2. The losses for the pump in the previous question are estimated to be 15% of the relative dynamic 
head at inlet to the impeller and 45% of the dynamic head at inlet to the diffuser. Calculate the 
head delivered and the efficiency of the pump. Answers: 4.12 m; 75.4% 
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Chapter 11 


Summary 


There are three really difficult ideas in this book: 


The Cascade and Meridional View. Essential this idea is how one can relate the machine 
geometry to the analysis model. In the author’s view this is one of the most difficult concepts 
in turbomachinery and is the most important idea in the book! This concept has been discussed 
at length throughout the book but with particular emphasis in Chapter 1. 


Velocity Triangles. These are essential for the calculation and understanding of machine per- 
formance and allow the graphical representation of the “key rule”: V = Ü + W, this concept 
is used throughout the book but is first explained in Chapter 2. 


Reaction. Very simply reaction is a choice between high work (low reaction) and high effi- 
ciency (high reaction). The geometry set by the blade angles determines velocities through the 
turbomachine and in Chapter 6 it was shown that this determines the reaction: the pressure drop 
between the stator and the rotor. 


A number of key equations have been developed and applied to turbomachines: 


e Continuity or conservation of mass. Equation 5.1: m = pAV 
e Conservation of angular momentum. Euler equation 5.12: w = w(Vogr2 — Vier1) 


e Conservation of energy or Rothalpy conservation across a blade row. Equation 5.17: I = 
ho — U Vo = constant 


Finally the concepts and equations have been combined to analyse the performance of a number 
of different machines: gas turbines, steam turbines, pumps and hydraulic turbines. This book gives 
an introduction to turbomachinery and (hopefully) empowers the reader to either understand more 
advanced texts or give an appreciation of the operation of these important devices. Further reading 
on turbomachinery may be found in Cumpsty (1997), Dixon (2005) or Cohen et al. (1996). 
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Appendix A 


Glossary of Turbomachinery Terms 


Adiabatic Process A process with no heat transfer 
Aerofoil An aerodynamic shape designed to produce a lift force with minimal drag (UK spelling) 
Airfoil An aerodynamic shape designed to produce a lift force with minimal drag (US spelling) 


Blade The name given to the part of the turbine, compressor or fan that guide the flow. Can be 
stationary or rotating. 


Blade Passage The space between two adjacent rotor or stator blades. Essentially this is the empty 
space that the fluid flows through. 


Blade Span The length of the blades in the radial direction 
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Bucket Another name for Blade 
Casing The higher radius where the case of the machine is located 
Chord The shortest distance between the leading edge and the trailing edge of an aerofoil 


Compressor Device that delivers energy to a fluid, usually used to describe a machine operating on 
a gas as the output is at a higher density to the inlet 


Diaphragm Another word for stator, usually used in steam turbine practise 


Enthalpy Formally: A thermodynamic property which is the sum of the internal energy of a sub- 
stance and the product of the pressure and the specific volume. In more practical terms it is a 
measure of the energy of a fluid at a given pressure and temperature 


Entropy A thermodynamic property that is constant for ideal process and increases for real processes 


Fan Devices that deliver energy to a gas. Usually restricted to a device at has a modest pressure rise 
over the device. 


Gas Turbine The name given to a complete machine including a compressor, combustor and turbine 
using air as a working fluid. Used to differentiate from a Steam Turbine which uses steam as a 
working fluid or as a generic name for the technology. 


Guide Vane Another name for stator blades, often used with hydraulic turbines 
Hub The lower radius at which the blades are attached 


Hydraulic Associated with liquids, for example a hydraulic turbine is one that uses liquid (usually 
water) as a working fluid 


Impeller Rotating element of a pump 


Isentropic Process A process with constant entropy - that is a perfect process that will have the 
highest possible efficiency 


Isothermal Process A process that takes place at constant temperature 
Isobaric Process A process that takes place at constant pressure 
Leading Edge The point on an aerofoil that is furthest upstream 
Nozzle Another name for a stator 


NPSH Net positive suction head - the margin between the saturated pressure vapour and the pressure 
at inlet 


Pump Another name for a device that delivers energy to a liquid. The term compressor isn’t appro- 
priate as liquids are almost incompressible. 


Rotor Rotating row of blades found after the stator in a turbine and before the stator for a compressor. 
Runner Name for the rotating element of a machine, often used in hydraulic turbines 


Scroll Casing A casing that resembles a sea shell for a radial turbomachine which aims to ensure an 
equal distribution of flow around the periphery by reducing the area of the device. 


SFEE The steady flow energy equation a fundamental equation of Thermodynamics. 


Stage A stator and a rotor combined, can be found either in a compressor or in a turbine. 
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Stagnation Conditions Conditions (pressure, temperature) that would occur if a fluid is brought to 
rest in an isentropic manner 


Static Conditions Conditions (pressure, temperature etc) which occur in a fluid which is moving 


Stator Stationary row of blades found before the rotor in a turbine and after the rotor for a compres- 
sor. 


Station (Or Analysis Station) A convenient point in a turbomachine where analysis can take place, 
such as the entry to a blade row or the exit from the machine as a whole. 


Steam Turbine A turbine using steam as a working fluid. 
Swirl Having a tangential component 
Tip The upper radius of rotor blades 


Tip Clearance The gap between the moving blades and the stationary casing in a turbomachine. For 
high performance aeroengines this is around 1% of the blade span. 


Tip Leakage The name given to the flow that goes through the tip clearance rather than following 
the path of the main flow through a turbomachine. 


Trailing Edge The point on an aerofoil which is furthest downstream 
Turbine Device for extracting energy from a fluid 

Vanes Another word for blades 

Vaneless Having no vanes or blades, used often to describe diffusers in pumps 
Volute Another word for scroll casing 


Wicket Gates Another word for guide vanes 
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Appendix B 


Picture Credits 


Figure 1.1 top left image is from Flickr user bcorreira. Licensed under the Creative Commons Attri- 
bution 2.0 Generic License accessed 30th November 2008. 
http://www. flickr.com/photos/bcorreira/2900837176/. 


Figure 3.1 left image is from Flickr user phault. Licensed under the Creative Commons Attribution 
2.0 Generic License accessed 2nd February 2009. 
http://www. flickr.com/photos/pjh/185488383 


All other images are from the Author’s collection. 
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